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ABSTRACT 


A new approach to numerically model two-dimensional confined flows 
is presented. The flow in the duct is partitioned into finite streams. 
The difference equations are then obtained by applying conservation prin 
ciples directly to the Individual streams. Also presented is a listing 
of a computer code based on this approach in FORTRAN IV language. The 
code computes two-dimensional compressible turbulent flows in ducts when 
the duct area along the flow is specified and the pressure gradient is 
unknown . 
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1 . INTRODUCTION 


Two-dimensional effects, such as the growth of momentum and thermal 

boundary layers, in flows through nozzles, channels, and diffusers have 

been traditionally studied by separating the flow into a core and a boundary 
1 

layer. The core flow is then calculated by one-dimensional potential flow 
equations and the boundary layer by momentum integral approach. Some 
flows cannot be accurately calculated with this approach. For example, 
in flows with large differences between the core and the wall temperatures 
substantial heat flow in or out of the core weakens the assumption of a 
radially uniform core. Again in calculating performance of diffusers 
with high inlet blockage factor, as in MHO diffusers, the flow cannot be 
meaningfully separated into a core and a boundary layer. 

Recently many approaches, based on finite difference methods, have been 
developed to calculate confined flows without the separation into core and 
boundary layer. These approaches can be broadly divided Into two classes. 

In one class the finite difference equations are obtained by differencing 
the flow differential equations either in their primitive form or in 
some transformed form. In the other class the finite difference equations 
are obtained directly from the conservation principles. Examples of 
these two classes, along with their interrelationship, are given in 
Fig. 1. 

This report deals with the approach of Ref. 7—shown bordered by 
heavy line in Fig. 1. The original paper describing this approach is 
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Included in Appendix B. In section 2 of this report are given the under- 
lying ideas of this approach, summary of the difference equations and 
a discussion of the solution technique. Based on this approach a computer 
code, called NEMCO (numerical modeling of confined flows), has been devel- 
oped. This code computes basically two-dimensional mass, momentum, and 
energy transfer. Thus it can handle flow which is either strictly two- 
dimensional— flow through nozzles, pipes, conical diffusers, discs, etc.— 
or is a two-dimensional approximation to the actual flow. An example of 
the latter is flow through a rectangular channel where the calculations 
are carried out along one pair of the opposing walls and the presence of 
the other pair is either completely ignored (case of large aspect rati^) 
or approximately taken into account by introducing some sort of effective 
width for the actual width between the ignored walls. 

In section 3 is given the program listing of the code NEMC0-U2-79 in 
FORTRAN IV language. The program listing given was derived, with some 
FORTRAN language modifications, from the various versions of this code we 
have used in our computations. To check the code listing included in this 
report we have run two sets of sample calculations— one for flow through a 
pipe and another for flow through a rectangular diffuser. The results of 
the sample calculations are also given in section 3. 


2. STREAHWISE COMPUTATION OF DUCT FLOWS 


In this approach finite difference equations are obtained by first 
part; cioning the flow in the duct into a finite number of streams and 
then applying conservation principles directly to the individual streams. 

To partition the flow into streams we draw a series of streamlines 1,2,...,N 
as shown in Fig. 2. The direction x is along the flow and z, measured from 
the duct wall » normal to it. Duct wall is chosen as streamline 1 and the 
centerline or the symmetry plane as streamline N. We define streamline 
j + l/2 as the streamline in the middle of j and j + 1. Streamline j - 1/2 
is defined similarly. We now define the various streams as follows: 

Stream 2; flow between the wall and the streamline 2 + 1/2, 

Stream j (N > J > 2) ; flow between streamlines j - T/2 and j + 1/2, 

Stream N; flow between the streamline N-1/2 and the centerline 
or the symmetry plane. 

Shape of a typical stream depends on the shape of the duct: in rectangular 
ducts the streams will be plane layers and in circular ducts the streams 
will be cylindrical shells. 

We introduce the following nomenclature: 

ip j’*' J mass flow rate between the streamline j and j + 1/2. 

ipj" 5 mass flow rate between the streamline j and j ^ 1/2. 

uj(x), Tj(x), pj(x), hj(x), Cj(x); velocity, temperature, density, 

enthalpy and constant pressure specific heat along the streamline j. 

Zj(x); distance of streamline j from the wall. 

Aj^; cross-sectional area of stream j. 
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^j+1/2’ between streams j and j + 1. 

^3+1/2^^^* *^3+1/2^^^’ effective viscosity and thermal conductivity 
(including turbulent contribution) along streamline j + 1/2. 

Ap; pressure difference between x and x + Ax. 


V 


- ft s 


V 


^j+1/2 


■i+1/2 


’"3+1/2 


where Py. is Prandtl's 
number. 


a,e,6; correction factors. 


Applying conservation principles to individual streams across x 
and x + Ax we obtain the following finite di fference . equati ons for velo- 
cities and enthalpies (see Ref. 7): 

Stream 2 

-U3 ^ ^ ’^2*'^ [4^2 “’^'2 ^2+1/2 ^1 ] 

= -ApA^^ + U3_-|“+ *^2lj^2 + aii’2 “ i^2+l/2 “ * (2-1) 

- hg 1~ + ^2+1/2) ^2 (1^2 ^’^'2 ■*■ ^2+1/2 1^1 ) 

("1“ ^2+1/2) ■*■ ^2(1^2 ^’^’2 ■ ?*^2+l/2 ” 1^1) 


= h 




(-^+^ 


2 , lti+ 


4 2*2+l/2 


) - ^ ( ^’*'2' 14+1/2) 


2 + 


+ ^(^+k 


2 V 4 2"2+l/2 


) ■*■ 2 ( 5 ^ 2 ^ ^'^2 ■ 2^2+ l /2) *' w2^1 


( 2 - 2 ) 
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In Eq. (2-2) represents enthalpy evlaluted at the wall tempera- 
ture. This program assumes that the wall temperature is given. When 
wall heat flux is specified see Ref. 7 for the needed modification. 


. Streams-^ 3,4->. . . ,N-1 


■“j+l * K+1/2^ M+1/2 * K-V 2 ] 




■ “j-1 


4 2 j-1/2 


] 


-ApAj^ + ^j+l["i“ 2 ^ 1 / 2 ]'^ " 2^j+l/2 ■ ’^ 3 - 1 / 2 ] 

^^j-l [” 4 “ ^j-1/2] * 

+ 

■ 2*^j+l/2J''' *'j [|*J ■'■ 1*0 ^j+1/2 l’'j-l/2] 


- •'j-l [-- ^ + K- 1 / 2 ] 

Vl [+ ^ ?”j+l/2] * 

" '> 0-1 [Y ^ h-1/2] " ^r- " Vj-1/2] 


+ "S _ 1 


h. 


2“j+l/2 F'j-1/2 

+ 


] 


[l*j" " l*o' C , ,2 + v: ] + ^ * vt,,,] 


4’^J ’'j+1/2 ‘ ''j-1/2 
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I ^ 

2 L 4 • ‘J+ 1 / 2 J 


2 . + 


+ V.., 


:] + ^[|co" + - V,. 


+ . 3 - 

^''j " ''j+1/2 " ''j-1/2 


4 ■ ‘j-1/2- 

] 


u • 1 
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+ V, 
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(2-3) 
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(2-4) 


stream N 


Difference equations for stream N are obtained from Eqs. (2-3) 
and (2-4) by letting j = N, and setting 

'*'N ^ \^^/2 ^ Vl/2 " “m+ 1 "" Vl “ 

^N+1 " ^N+1 ^ ^N+1/2 " *^N-M/2 " ° 

In the equations superscripts + indicates the quantities at x + Ax, thus 

Uj is the velocity along j at x and Uj along j at x + Ax. In the listing 

given in section 3 same symbol is used for a given variable, like uj, 

at x and x + Ax. Program is written such that from the point calculations 

of new quantities are started there is no further need for old values 

at X. Equations (2~1) to (2-4) are supplemented with formulae given 

in Appendices A and B of Ref. 7. These formulae calculate A^, A®, li, 

+ 

and first at x calculate tp- and if . given p., u., and z- and then later 

J J w vl li 

+ + - + + 
at x+Ax calculate z. from ip. , i|>. , p 4 , and u. . 

V «J J J J 

The technique to solve these difference equations depends on the 
specified design variables. The program given in section 3 assumes that 
the variation of duct area along the flow is specified and the pressure 
gradient is unknown. For this case the difference equations are solved 
as follows. Equations (2-1) and (2-3) (and similarly Eqs. (2-2) and 
(2-4)) are looked at as a system of tri-diagrnal linear algebraic equations. 
This, of course, is true only if V.^ and Ap are known. To start with 
we set V. equal to V. and solve for u. using standard methods to solve 

J J J 

“t* 4* 

tri -diagonal equations. We then calculate V. from the new u. and solve 

<J J 


for u. again using the new V. . This iteration process is continued 

J J 

until new u-^differs from the previously calculated by less than 

J J 

some preassigned value. In the early computations u«» among all u., 

C J 

was found to change most rapidly from one iteration to another. Thus 

it has been chosen for all our computations to check for convergence. 

This iteration loop is called the velocity iteration loop. Solution for 

h. is carried on along with the velocity iteration calculating H. from 
•J J 

the latest V. and using the latest available u. for kinetic energy 

V J 

terms. 

The pressure drop Ap is also not known in advance. We start the 
calculations with some estimated Ap and carry out the velocity iteration. 
When the velocity (U2^) convergence is achieved we check the calculated 
duct area at x+ax against the specified duct area. If the difference 
is larger than some preassigned value we go back, adjust Ap, and start 
the velocity iteration all over agian. Iteration on Ap is called pres- 
usre iteration. Note the pressure iteration loop is placed over and above 
the velocity iteration loop. Initial estimate of Ap at the beginning of 
a new integration step is obtained by projecting the previous three 
values using least square method as discussed in Appendix A. Correction 
to Ap from one pressure iteration to the next is based on the difference 
between the calculated and the specified duct area. 


Mr 
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3. PROGRAM LISTING 

3.1 Introducticn. 

The program consists of a MAIN PROGRAM and nine subroutines; 

INPUTl 
START 

SIPSIM 
ZDZ 

XSAREA 
VSCSTY 
TDSOLV 
THICK 

HPTRHO 

The main program and five of the subroutines are common to all duct shapes. 
The four subroutines INPUTl, SIPSIM, ZDZ, and XSAREA depend on the shape 
of the duct. Subsections 3.2 through 3.6 contain computer listings in 
FORTRAN IV as follows: 

Subsection 3.2; Main program and subroutines common to all duct 
shapes, 

Subsection 3.3; The four subroutines for circular ducts. 

Subsection 3.4; Sample computations— pipe flow. 

Subsection 3.5; The four subroutines for rectangular ducts. 



(Provides geometric data for the duct), 

(Sets up grid normal to the flow and assigns starting velocity, 
temperature, density and enthalpy profiles), 

(Calculates and , the mass flow rates through upper and 
lower parts of streams), 

(Calculates distance of streamlines from the wall), 

(Calculates cross-sectional and surface areas of the streams), 

(Calculates laminar and turbulent viscosities), 

(Solves tri -diagonal set of algebraic equations), 

(Calculates various thicknesses associated with the boundary 
layer) , 

(Calculates temperature and density given pressure and enthalpy). 


II 


Subsection 3,6; Sample computations—rectangular diffuser flow. 

An computer listings are proceeded by definitions of the FORTRAN names 
used in the program. 

DATA INPUT : 

Data input to program is through MAIN program and subroutine INPUTl. 
Input statements have been placed at the very beginning of the programs. 
Data given through the various programs is as follows: 


Main Program : 

a) UCNTR, ZBL, TCNTR, TWALL, P, CP, GASK, GASR, PRNDL, DPDX— 
dynamic and thermodynamic flow data. 

b) XEND— length of the duct. 

c) DPDTL, DXFAC, NJ, MAXIT, NPRNT—pressure gradient convergence 
tolerance, integration step size, number of grid points nonnal 
to the flow, maximum number of iterations allowed, and output 
frequency. 

d) VZERO, TZERO, EXPV—parameters needed in subroutine VSCSTY to 
calculate laminar viscosity. 

e) EXPU, EXPH- -parameters needed in subroutine START to calculate 
initial boundary layer profiles. 


INPUTl Subroutine : 

Duct shape, parameters differ for different duct shapes. 


PROGRAM OUTPUT 


Output of program is of three kinds: 

a) Distribution of velocity, temperature, density, and enthalpy 
profiles normal to the flow at a given value of x. Frequency 
of this output is governed by NPRNT. 
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b) Distribution of pressure* shear stress, heat transfer to the 
wall, and various boundary layer parameters along x. These 
quantities are calculated at the end of each integration step 
and stored in vector forms. The results are printed after the 
calculations reach the end of the duct. 

c) Many subroutines have output statements to print results of 
intermediate calculations. Frequency of this output is gov- 
erned through the subscript appended to the subroutine nsrnie. 
This kind of output is used for diagnostic purposes. 


3.2 Listing main program and common subroutines. 


MAIN Program ; 

U(J) = Velocity along streamline J 
TfJ) = Temperature along streamline J 
H(J) = Enthalpy along streamline J 
RHO(J) = Density along streamline J 
USH(J) = (1/2)U2{J) 

Z(J) = Distance from the wall of streamline J 
DZ{J) = Z(J+1) -Z{J) 

NJ = Number of grid points normal to the flow 
NO! = Nd-1 
NJ2 == NJ-2 
m = NO-M 

NBL = Number of grid points assigned to the boundary layer at the 
inlet (NBL<NJ) 

Note: In this program NBL is calculated internally in subroutine START 
DX = Size of Integration step along the flow 
SlP(d) = 

SIM(J) == 

FLORAT = Total mass flow rate through the duct 

UCNTR = Velocity at the .centerline or the symmetry plane 

TCNTR = Temperature at the centerline or the symmetry plane 

TWALL = Wall temperature 

P = Pressure 

CPINV = Inverse of Cp, the constant pressure specific heat 

RINV = Inverse of R» the gas constant 

fiASK = Ratio of the specific heats 

ZBL = Boundary layer thickness at the inlet 

AX(J) = Cross-sectional area of stream 0 

AS(jj = Surface area between streams J and J-1 

ALPHA, BETA, DELTA, GAMMA are the correction factors ct, 0 , 6 ,Y-of Ref. 7. 

RE Reynolds number along streamline 2 

VISC(J) = Viscosity along, streamline J-1/2 

DELI - Displacement thickness of the boundary layer, 

DEL2 = Momentum thickness of the boundary layer, 62 
H12 = Shape factor DEL1/DEL2 
DELO = Boundary layer thickness 

C0MM0N/C0M8*, this common statement contains variable names pertaining to 
the duct shape. These names are defined later in INPUTl. 

VF(J),HT(J) defined later by the program listing. 

UU(J) = Intermediate velocity variable 


The tr1 -diagonal form of the momentum and energy difference equa- 
tion is 


AA{J) »BB(J)sCC{J) are used for aj, bj, and cj, first for momentum and then 
for energy equation; DD(J) is dj for momentum equation and TT(J) dj for 
the energy equation; DU(J) is part of DD(J) and DT(d) part of TT(S). 

Note: Subscript I, in what follows, denotes the integration step along x. 

EXVEC(I) = Distance along the flow, "x“ 

PVEC(I) = Pressure along x, 

SHRVEC(I) = Shear stress at wall along x, 

HTVEC(I) = Heat flux along x, 

DLOVEC(I) = Boundary layer thickness along x, 

DLIVEC(I) = Displacement thickness along x, 

DL2VEC(I) = Momentum thickness along x, 

SFVEC{I) = Shape factor along x. 

CP = Constant pressure specific heat 
6ASR = Gas constant, R 
PRNDL = Prandtle number 

DPDX = Pressure gradient guess at the entrance 

Note: The program requires as input a value of pressure gradient at the 

duct entrance for its first iteration. We suggest a value calculated 
approximately for one-dimensional flow. 

XEND = Length of the duct 

DPDTL = Tolerance within which dp/dx is to be calculated 
DXFAC = Specifies the size of integration step, DX = HEIGHT/DXFAC 
MAXIT = Maximum number of iterations allowed, if this number is exceeded 
the program will stop and write a message 
NPRNT = Governs how often along x U(J),T{J), etc., are printed 

VZERO,TZERO,EXPV: See subroutine VSCSTY 
EXPUaE XPH; See subroutine START 

XAMP,DXF,DX; for the first integration step DX Ts set e^ual to D^F/^56, 
where DXF is calculated from DXF=HEIGHT/DXFAC. DX is multiplied 
by XAMP at the end of each integration step when DX reaches 
DXF XAMP is set equal to 1 . 

ISTEP = Integration step counter 
ITRN = Iteration step counter 
NPI = Counter for NPRNT 
SNDSPD - Sound speed at the inlet 
PTINL = Total pressure at the inlet 

OLDA - Cross-sectional area of the duct at the previous integration step 
DPDAMP, ZTLRN, UTLRN; See Appendix A 

Symbols appearing in the rest of the MAIN program are either already 
defined or defined by FORTRAN statements appearing in this part of the 
program. 


START Subroutine: 


This subroutine sets up grid normal to the flow and specifies velo- 
city, temperature, and density profiles at the beginning of the compu- 
tations. The listing given assumes that the flow at the inlet has a 
boundary layer of thickness ZBL and that within the boundary layer velo- 
city and enthalpy vary as 

Jilil = and 

H(J1 -HWALL _ 

HCNTR - HWALL “ V ZBL / 

where HWALL and HCNTR refer to enthalpies calculated at wall and center- 
line temperatures, and EXPU and EXPH are input parameters. Outside the 
boundary layer H and U are assumed to be uniform. 

To set up grid that will compute the flow with maximum possible 
accuracy is difficult. Algorithm given in this subroutine is guided by 
the following considerations: 

1. Grid points are to have greater density near the wall. 

The grid spacing to increase as Z to the some power, where 
Z is the distance from wall. 

2. Grid point Z(2) (Z(l) is wall) is to be located just outside 
the laminar sublayer and the transition zone — the formula used 
to calculate turbulent viscosity for points Z{J), 0>2, is good 
only for the turbulent region. Also point Z(2) must not be 

too far away from the wall or we will lose accuracy in calculat- 
ing flow near the wall. 
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3. Grid point Z(3) is to be selected so that Z(3)~Z(2) is not 
less than Z(2)-Z(l). Otherwise, as noted in Ref. 8, for 
some flow conditions one may encounter numerical instability. 

In the subroutine Reynold's number along Z(2) is required to be around 
1200. Thus U(2) and Z(2) are calculated from the two simultaneous equations. 


UCNTR \ ZBL 


/EXPU 


and 


U(2) Z(2)/VISKIN = 1200 


. where VISKIN is kinematic viscosity. 

Grid location of Z(J), J&2, are next calculated using a variation 
of formula given in Ref. 8. We put 

/.nEXPG . 

Z(J) = Z(2)+ (HEIGHT -Z(2))(^j 
Grid exponent EXPG Is calculated from the requirement 


Z(3)-Z(2) = Z(2)-Z(l) 


VSCSTY Subroutine; 


This subroutine calculates viscosity (laminar + turbulent) coeffi- 
cients along the stream interfaces. Symbol VISC(O) represents the vis- 
cosity coefficient along streamline J-1/2. The laminar contribution to 
VISC(J) is calculated from the formula 


VISC(J) = l[vZER0 + VZERO 



where VZERO 3 EXPV, and TZERO are constants and are required input to 
the subroutine. The turbulent contribution to VISC(J) is calculated 
from 


VISC(J) = r a(J) + EL(J-1)f rRHQ(J)-fRH0(J-1) 


DZ{J) 


The mixing length, EL, Is calculated from the relation for the fully 
developed turbulent pipe flow [Eq. (20.18), Ref. 1]. 

EL(J) = ZOUT [o.l4-0.08 p 0.06 ^^^^^) 

The quantity 1 -Z(J)/Z0UT represents the distance from the centerline 
normalized by half the channel width. 

This subroutine also calculates the correction factors y, a, fi, 
and 6. The formulae used in the subroutine to calculate y* and 5 
were obtained by curve fitting the graphs shown in Ref. 7. Factor & 


was calculated from ct using the relation given in Ref. 8. 
is obtained using integrals defining tx and p, and assuming 


(This relation 

^ ^ JL ). 

h " hg ug 


The formula to calculate FF from Re was obtained by curve fitting the 
relation shown in Fig. B-1 of Ref. 7 in the vicinity of RE- 1000. 
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TDSQLV Subroutine : 

This subroutine solves the tri -diagonal set of equations 







where j = 2,3,...,N. Method used is the well-known special case of 
Gauss' reduction formula (see for example Ref. 2). The boundary condl 
tion at the wall is already built into the difference equations and 
appears as 

C2 = 0. 


The boundary condition at the center line used in the present work 
is that of symmetry and appears as 


Ah = 0. 


THICK Subroutine; 


This subroutine calculates boundary layer thickness, fig* displace- 
ment thickness, 6*], momentum thickness, 52 » shape factor H.jg, 

Boundary layer thickness is defined as the distance from the wall 
where the velocity equals 0.995 times center line velocity. The integrals 


CO 



were approximated by the sums 


(' • i J 

(i-^]h 

V P«Uo,/ 

V P<»Uo<, / 

j=3 


j 


and 



\ PcrjUgo 


^TP2U2 

pcoUto^ 





PoqUco V Ura ) Poo^co \ 

u„ i 

J-3 

- 
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HPTRHO Subroutine ; 

This subroutine calculates temperature and density given pressure 
and enthalpy. In its present form it uses formulae for perfect gas 
with constant specific heats. 


i FUKTRAN 
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IV 360t^FU-479 3-8 MAINPGM DATE 08/20/79 TIME 13. 

C 

!(««:'«=***** :4t* sr si- *<:*#**** 

C 

CQMMQM/C0M1/UI50) »T tSOJ jHI 50} ,RH0C50J ,USHI50) tZ { 5 0} t D2 t 50 ) 
COMMON/COM2/NJ»NJlTNJ2rNJPf NBLfOX 
C0MHQN/CDM3/SIPI50} tS IH ( 50 J » FLORA! 
CGMHaN/C0M4/UCNTRTTCNTR,TWALLTP,CPINV,RlNV7GASKt28L 
C0MM0N/CQM5/AXI50 1 , AS 150} 

CUMM0N/C0M6/ALPHA,BETAt DELTAtGAMMA, RE ,VISC{ 50 ) 

COMHON/COM7/DEL1 tDEL2TH12,0EL0 

CaMMON/CQMR/HEIGHTr WIDTHtHSLDPEjWSLCPE, AVWlDHjEFW IDHt CLOA ,APF AC 
CQMM0N/CDM9/V2ER0,EXPV» TZERQ, EXPLf , E XPH 
DIMENSION VFt 50J ,HT{50) tUUl 50) 

DIMENSION AA( 50 J »BB{ 50) , CC £ 50 } , OD ( 50 ) tDU ( 50 ) r DT ( 5 0 } » T T ( 50 ) 

DIMENSION EXV£C£200),PVEC£200J»PRVECI20J) tSHR VEU 200) , 

IHT VEC £200) r DLOVEC £ 20 0 J t DL1VEC£200) i 0L2VEC 1 200} , SFVECf 200 ) 

C 

C'S^^DA TA I NPUT **>};>{;***:{: ip + ^Issit J):s^ ^:{c ft# *!{■*« #ft##ft# #1)* Aft*######### 

C**DVNAMIC AND THERMODYNAMIC ##### 

UCNTR=SM 

ZBL=«aM. 


CP = 

GASK=1 
GASR-I 
PRNDL=| 

OPDX = i 

C*+GEUHETRIC DUCT 
CALL INPUT! 

xemd=ms 

C=t'*CUNVEKGENCE CR I T EK I A f GP I D SIZE, AND OTHER o AT A” *>>' ^« ♦**>*'** ^ 

DPDTL^I 
DXFAC- 
MAXIT=i« 

NJ = WI 
NPRNT=»8 

C*=f=OATA FOR VSCSTY AND START SUBRCUTI NES****«’'**«***«*^!*'!‘^'^^* = 

VZtRO^I 
TZER0=1 

£XPV=aii 

fcXPU=lW 

EXPH=W 

C 

C>*'*INITIALIZATI ON S*********^^ ###***#* #*#* #>?:## ****** **#» *****»■ *ft*ft - ■■■ ■ 

XAMP=2.0 
IS TEP=0 


ITRN=Q 

NPI=NPRNT“1 

NJP=NJTl 

NJI=NJ-1 

NJ2=NJ-2 

SIP £NJ)-0. 


0 
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U fURTRAM IV 360N-F0-479 3-8 


MAINPG^i 


OATF Oa/20/79 Tl^f 13 


( ■ 


i 


VF£NJP)=0. 

HTINJPJ=0. 

JINJP)=0. 

HINJPJ=0. 

USH(NJPi=0. 

CC(2J=0. 

OXf=HEIGHr/DXFAC 

DX=DXF/256. 

ex=-DX 


DP = DPDX=«'DX 

PRLIMV=1./PR(\|D1- 

CPINV=U/CP 

RIINJV=1./GASR 

GASKIN=1./GASK 

HWALL-CPiPTWALL 

SNJSPD-fGASK«GASR=PrCNTR )=P*0.5 

PTINL=P’!'( l.i-.5*IGASK“l. )’!'{ UCNTR /SNOSPD ) 2 )’!'*( GAS K/ t CASK- L . ) I 

PIf\lL=P 

DPPTI=PriNL-PINL 


DPi=0, 


DP2=0. 


0X1=0. 


OLOA=0.O 
CALL START 
WRlT£f6 ,30} 

30 FCiRMATl IHl, « VELOCITY U , TEMPERATURE T, DENSITY RHC), AfJO ENTHAPY it 
iDISTRIElUTIQNS AS FUNCTION OF Z.*,//) 

C 

C^^AX lAL CALCULAT 
1100 ISTEP=ISTEPF1 
EX=EX+DX 
CALL VSCSTYIO} 

CALL THICK 
NPI = NPI-f-l 
EXVECnSTEP J=EX 
PVEC I ISTEPJ =P 

PRVECI ISTEP}=iP-PINL)/DPPTI 
SHRVeCUSTEP}=VISC( 2)=PU(2} 

HTVECnSTEPl = VlSCI2 )=J=PRLINV’PIH( 2i-HWALL} 

DLOVECI ISTEP} =DELO ■ 

DLlVECI ISTEP}=DEL1 
DL2VECI [SrEPJ=OEL? 

SFVEC(ISTEP)=hL2 

IF! NPI.LT.NPRNTJ GO TO 1200 

NPI=0 

WkITE(6 ,46} EX ,RF 

46 FORMAT! **^^*OIST- ALONG X= ' , F7. 4 , • ^*NE AR WALL R E= • , 1 PF 1 0 .4 ) 
HRITEI6,40) (Z{ J} ,J = 2,NJ) 

WRITEl6,41i { Ut J) ,J=2 ,NJ} 

WKITE(6r42) (T(J},J=2,IJ} 

WRITE(6,43J { RHO( J3 , J=2 ,NJ) 

WRI TE (6 ,44) (HtJJ , J=2,NJJ 

40 FORHATI/tlX, *Z‘ ,aX , IP lOEi 1.4/ { 1 OX ,iPiO£il .4 J ) 

41 FORMAT!/, IX, «U* t8X , IP 1 OEl 1 . 4/ ( 1 OX , 1 P 1 OE 11 .4 J ) 
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I ; 

1 I 

s 

IjS FORTRAN IV 360N-FO-479 3-8 MAINPGM DATE 08/20/79 TIME 13 

1 . : 

;■ 42 FOKMATf/,lXT*T' t3X t 1? lOEIl ,4/ ( lOX , i PlOE II -4 ) 3 

; 43 FORMAT(/t* RHO' ,6X » IP iOEl 1 .4/ ( lOX, IPlOE ll .4 ) } 

44 F0RMAT(/t 1X, »H» ,8X , IP lOEll, 4/ ( 10X,1 PlOE I I .4 ) 3 

! IF (EX.GT.XEND] GO TO 2L0Q 

C 

RESETTING DYNAM IC tTHERMOf AND GEOM, VARIABLES FOR NEXT I T ER AT I 
1200 ITRN=1 
Xl=OXl 
X2=DX+0XL 
SUMX=X1+X2 
SUMXX~X1=«‘XI+X2=S=X2 
i SUMP=0P+DP1+DP2 

u SUMXP=XJ^QPH-X2#0P 

f: DE M=3.*SLfMXX-SUMX^SUHX 

f. OPSLOP={3.^SUMXP-SUMPi'SUMX3/DEN 

' DP2ERQ=ISUMP4SUMXX-SUMX*SUMXP3/DEN 

j; DXi = DX 

DP2=DP1 

t OPL=DP 

1F( ISTEP-LE.2) GO TO 1205 
! IMDX-GE.QXF3 XAMP=UO 

f DX=DX’>'XAMP 

IFI ISTEP.lt. 43 GO TO 1205 
OP-DPSLGP*( X2+DX) +DPZERO 
1205 PS1ART=P 

UPJAMP=a.5=J=EXPt-0.1l5=S=EX/HE IGHT 1 
U0LD=U{2) 

; CALL SIPSIMtO) 

C CALL XSAREA(O) 

ZFAC=FLORAT*U(NJ)/I OLDA-FLORAT^UiNJ )=pGASKIN/P 3 
ADP=ABS10PL3 

r ZTLRN=ADP^=DP0TL4HEIGHT/ZFAC/ARFAC 

UT LRN=ADP>?‘0PDTL*0LDA/ FLORA r 

i: c 

C=!=*-=«= EVALUATION DF K.H.S. OF MQM. AND ENERGY EQS. ’J-" ** « 

VFI 23 = .5=<=A5i2)*VISC(2) 

^ HTi23=VFl2)*PRLINV 

DO 1106 J=3,NJ 

VF I J 3 = .5>s=AS ( JI^'VISC U 3/DZI J 3 
t 1106 HTI J3=VF{ JJ^^PRLINV 

DUI 2 3=0 (3 3^ (S IP t 2 3+VFI 3 3 3 i-U { 2 3 * ( 3 . S I P I 2 ) + ALPHA*' S IM ( 2 1 - VF (3 ) - VF { 2 ) 

1) 

DD(2)=DU(23-DP’«'AX(23 
i: DO L108 J = 3,NJ 

i; L3U( J) = UIJi-l )^(SIP (J 3+VFU + l) )+UIJ3>!=(3.=(‘SIP( J) +3.*SIMU1 -VF( J^-l) 

1-VFC J3 )+U(J-l3»ISIM( J)i-VF( J3 ) 

DDi J)=OUU)-DP*AX{J} 

1108 CONTINUE 

DT<2 )=H(3 3*= (S IP [2 3 + HT(3 3 )+Ht 2 23 +BETA*SIM( 21 -HT{ 33 -HT( 2) ) 

1 + USH(3)=J=(SIPI23+VF[ 3 3 3 ^US H { 2 3 «S. IP ( 2 I +OFL T A=!=S I M{ 2 1 -V F t 3) )+HW ALL 

2’fHT[2) 

DO 1.112 J=3,NJ 

; DTI J3=Ht J-(-l)'4'(SIP { J3 + HT ( J + 13 3 *HU3=i'{3.*SlP{ JI+3.=i=SIM( Jl-HTiJ-i 13 

i i-HTC J3 H-HIJ-D^tSIMIJli-HTt J3 3^-USH[ J-M)=«'(SIP { J3+VFI J + 1 3 3 +USH( J3* ( 


sss.«- 
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FORTRAN IV 360N*-FQ-479 3-8 


MAINPGM 


DATE 08/20/79 


23-=S=SIP( J]+3.=i=SIM[ J}-VF( J+IJ-VFC Ji ) +US H ( J~1 } f S IN ( J ) +VF tJ ) ) 
1112 CONTINUE 


C^e^ITERATIVE CALCULATIONS FOR NEW U »H rZ 

C**CALCUL AT IONS OF NEW VELOCI TI £5*’^^*******+'^*^=?'=^*********** 

GO TO 95 
90 CONTINUE 

VF(2I = ASI 2)’f'VISC(2J*0.5 
KT12}=VF(2)*PRLINV 
DO 99 J=3,NJ 

VF( J3=AS(JJ *VISCI J)/OZ( Ji^=0.5 
99 HT£ J)=VFIJ)^PRLINV 
95 AA(2)=-SIP(2J+VF(3} 

B8(2J=3.^SIP(2J+ALPHA^S 1M{ 2 ) +VF { 3 ) +VF (2J 
DO 100 J=3tNJ 

AA £ J » =-S I P ( J J «-V F ( J + n IS 

B0£ J)=3.*SIP£ J)+3.*SIM( JJ+VF£ J<-L)+VF£J) 

CC£ J)=-SIM£ J)+VF£J) 

100 CONTINUE 

CALL TDSOLV ( AA ,BB »CC »3DtUU ,0 ) 

IF HTRN.GT.U go to 109 
DO 108 J=2,NJ 
100 UtJJ=UU(J> 

109 DO 200 J = 2,NsJ 
200 UIJ >=.5=<={U( J) +UU( J) 3 
C 

Cv+CALCULAT ION OF NEW EnTHALP 
DO 111 J=2fNJ 
111 USH£ J 3 = .5=»^U( J)*U£ J3 
AA£23=-SIP(23+HT£3) 

BB( 2)=3.*SIP£ 23+8ETA*SIiM(2M-HT( 33 +H r£2) 

TT £2 3=DT{2)+USH(3)*£-SIP(2) +VFt 3) )-USH{2) ”£ 3.’>=SIP ( 2 3 +.)t I T S I -U 2 > 
U-VF(3) ) *HWALL-HT£2) 

DO 115 J=3»NJ 

AA{ J3=-SIP£ Jli-HTf J^-U 

Bd£ J)~3.=i=SIPf J)+3.=S'SrMI J)+HT£ J + L) +HTIJ3 
CC( JJ=-SIMf J) +HT£ J3 

TT( J)=OT£ J3+USH(J4-1 3*£-SIP£J3+VFfJ + U 3-U5HU3 -( 3. «S I P t J 3 + 3. =? S IM ( J ) 
l+VF ( J + 1 ) + VF( JH +USH( J-U«£-SlMt J) frVF£J3 3 
115 CONTINUE 

CALL TDSOLVIAAtBB rCC,TTTHt03 
C 

P=PSTART+DP 
CALL HPTRHO 
CALL ZD2 

C^^CUVERGENCE CHECKS 4- 

ZCGK=Z£NJ) -HEIGHT 
80 ITRN=ITRN+1 

UCHEK=ABS(U(2)-UaL03 
ZCHEK=ABS£ ZCOR) 

IF £ UCHEK.GT.UTLRN3G0 TO 82 
IF£ZCHEK.LT.ZTLRN3G0 TO UOO 
IF ( ITRN.GT.MAXIT) GG TO 2000 
DPCOR=-ZCCR*ARFAC/Z £NJ3 =^=ZFAC 
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- ;S FORTRAN IV 360N-FQ-479 3-8 MAINPGM DATE 08/20/79 TI^IF 13 

310 DPCOR^DPCQR^i^OPDAHP 
330 DP=DP+DPCOR 

I ba 340 J=2,NJ 

340 DD{J}-OU{J)-DP=^AX(Ji 
r 82 CONTINUE 

V UQLD-UI2) 

; ■ CALL VSCSTYtOl 

i GO TO 90 

: c • 

. > 2000 WRITE(6,2D0U 

1; 2001 FDRMAT{/TlX,»4*VELaCITY-PRESSURE ITERATIONS EXCEEDED THF LIMIT SET 

r t' 1 by maxi t*******4***^* ) 

2100 WRITE(6f21011 

2101 FORHATt IH1»9X,*EX* ,9 Xt' PRESS* t9Xt 'P-RECV* t6X» ‘WALL SHEAR* 

U 1,3X»*HEAT FLX',5Xt* fl-L. THK ' , 5X, ' 01 SP THX ' , 3X , • MOM Tt^K * , 6 X, * S MAPS 

V' 2EAC*,/) 

WRITE(6,2105J ( EXVECIK) ,PVEC(K) fPRVECtK} ,SHRVECtKJ , 

IHTVECtK) ,OLOVEUKI,DLlVECtK3,OL2VEC(K ItSEVEC(K) , K=1,ISTEP) 

2105 EaRMATtlPEl6.3,LP8El3.3) 
i STOP 

ENO 






START 
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FORTRAN IV 360N-Fa^A79 3-8 


DATE 08/31/79 T 18. 


SUBROUTINE START 

C**SUBRQUTIN£ SETS INITIAL PROFILES AND NORMAL TO THE FLOW gR 

COMMON/CQMl/UfSO) ,T(5Q1 ,H ( 50) tRHO ( 50 1 » USH { 50} ,Zt50) rDZ(50] 
CaHM0N/C0W2/NJfNJI,NJ2r NJP tNBLiOX 

C0MMQN/C0M4/UCNTR,TCNTR»TWALL,P,CPINV,RINVtGASK,ZBL 
COMHON/COMS/HEIGHT, WIDTH, HSLOPEfWSLOPE, AVWIDH tEFW IDH tOLDA ,ARFAC 
CUMHQN/,C0M9/VZERa,EXPV, T2ERC, EXPU,EXPH 
C 

VISKIN=VZERa*(TCNTR/TZERCJ>!==i=EXPV/(P=!«RINV/TCNTR) 

FNJ=^FLGAT(NJ) 

HCNTR=TCNTR/CPINV 
HWALL=TWALL/CPINV 
D£LH=HCNTR-HWALL 
U2Z2=1200.=«*V1SKIN 
UU) = 0. 

z(u=o. 

DZ( 1)=0- 
H/U = HWALL 

U{2)={UCNTR=f'^=EXPU*U2Z2/ZBL)*>f={ I ./t l. + eXPUn 
Z(2)=U2Z2/U{2) 

DZI2)=Z(2J-ZU1 

HI2J = HWALL+DELH'{=( Z(2)/ZBL) **a./EXPH) 

EXPG=AL0G{HEIGHT/Z(2J-1 .)/ALQG(FNJ/3. ) 

DO 20 J=3,NJ 
FJ=FLOAT{ J) 

ZI J) = Z(2J -f-IHElGHT-ZIZJI^IFJ/FNJ I^^^^EXPG 
20 DZIJ)=Z{ J)-ZI J-IJ 
DO 23 J=3,NJ 

U( J) = UCNTR«(Z(J)/ZBL) I ,/EXPU J 

HI J)-HWALL+DELH=<^f Z( JJ /ZBLj=4=«{ U/EXPH) 

IF(UI Jl.GE.UCNTR) GO TO 25 
23 CONTINUE 
GO TO 27 

25 DO 26 K=J,NJ 
U(K)=UCNTR 

26 HIK)=HCNTR 

27 CONTINUE 
DZINJ+1)=0. 

PRINV=P*RINV 
DO 223 J=2,NJ 
TIJ)=H( J)*CPINV 

223 RHOt JJ=PRINV/TIJJ 
Tt D^TWALL 
RHUCl)=P*RINV/Ttl) 

DO 31 J=l,NJ 
31 USH( JJ=U(J}*U( J J*0.5 
RETURN 
END 
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IS FORTRAN IV 360N-F0-479 3-b VSCSTY DATE 08/20/79 


SUBROUTINE VSCSTYINPJ • 

CDNMQN/CONl/UtbO} ,T(50} tHiSO) tRHD(50) ,USH(50) tZ{ 50) tOZfSO ) 
C0MM0N/C0M2/NJf NJlfNJ2t NJPf NBLj DX 
C0MM0N/C0M6/ALPHA,BETA, DELTA » GAMMA, RE » V ISC ( 50) 
COMMQN/CQM9/VZEROtEXPVTTZERO,EXPU,£XPH 
DIiMENSrUN EL£50) 

C=S==«‘NIXING LENGTH , EL, CALCULATIGNS=5'***^^«**=«==«'^«’«'« 

ZOUT=Zf NJl 
E L f U =0 . 


DO 210 J=2,NJ 
FAC=(1.0-Z tJ)/Z0UT)**2 

EL { J ) = ZDUT«£ 0- L4-0.08=i'FAC-0.06*FAC*F AC) 

210 CONTINUE 

C*>5=VI SCQSITY ^LAMrNAR^-TUR3ULENT ) CALCULAT IONS*={==!==^'=i' 

VI set U=VZERO*tT( 1)/TZER0)=^*EXPV 
VI SOLD =VZERU*{Tt2)/TZ£R0)**EXPV 
RE=U(2)=f=RHO(2J*Z{2)/VISOLD 
D0215 J=3,NJ 

VISNEW =VZERO*tT( J)/TZER0J*=5'EXPV 
VlSL=0.5=S'fVISUL0i-VlSNEW) 

VISOLO=VISNEW 

ELrS=( tELt J J-i-ELt J-1 J ) *0.5)=s=«2 
VIST=ELTS*tKHO( J)+RHOtJ-l) j’f'O.S 
DUJ=ABS(U( J)-Ut J“l) ) 

VIST=VrST«DUJ/OZlJ) 

Vise tJ)=VlSL+VIST 
215 CONTINUE . 

FF=0. 0290 106/RE*« 0.2827 065 
VlSC{2)=FF*HHat2)*U£2 ) 

G=ALOG£RE) 

G2— G 
G3=G2=i=G 

GAMMA=0.6995594-0.1 7352 47*G+0. 043 43642«G2“0. 002524061*^3 
ALPHA=0.83l 1825-0,2 5511 32*G-«-0.0 5791 319*G2-0 .00 322 67 84*0 3 
0ELTA=0.79a7116-0,l 1948 7 7*G + 0.0 J0 75 32 7*G2-0.00180 62A1‘^G3 


i 

t 


102 

103 

999 


BbTA^tTt 1 )4-ALPHA* (T {2 )-T £ I J ) ) /T 1 2 ) 

IF £NP.NE .1 ) GO TO 999 

WRITE (6,102) (EL £Jl,J=l,NJ ) 

WRITE (6,103) (Vise £J),J=1,NJ ) 

FURMATf/r* £L',7X, 10£ 1 L . 3/ U OX , 1 OE 1 1 , 3 ) ) 

FORMAT(/,» VISC',5X, lOEl 1 . 3/ ( LOX , 1 OE L 1 . 3 ) ) 

CONTINUE 

RETURN 

END 


TIME 13. 
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•MS ^URTRAN IV 360N-F0-479 3-8 TOSOLV DATE 03/20/79 TIME 

} : 


i : 

I 








c 


SUBROUTINE TDSOLV I A »B rC , D ,X , NP ) 

.C0MM0N/C0M2/NJ»NJl,NJ2f ^^JP,NBLt DX 

DIMENSION AI50J tBiSO) ,C5 50)tD{50J ,X(50J ,E(5a) ,FI50) 
ROUTINE TO SOLVE TriE TRIDIAGONAL EQUATIONS 
E(2}=AI2)/S{2i 
F(2J=D(2J /Q(2i 
DO 500 J=3,NJ1 
DEM=B{J)-Ci JJ*£{J-l ) 

EIJ)=A{ J)/DEN 

500 FIJ)={D{JJ + C( J)^F{ J-UJ /DEN 

5 05 XlNJl )=IECNJ1 J^DINJ )FB( NJ)*F{NJ1J )/IBINJ 1-E(NJI)=}=C(NJ } ) 
DO 510 I=2 tNJ2 
J:=NJ-I 

510 XUJ=E(J)^X( J + U+FIJ) 

XiNJ5 = (0{NJ J+C<NJJi^XI NJin/BlNJ ) 

IF INP.NE.n GU TO 999 
WRITEl6,ia2I I AIJI t J=2»NJ1J 
. ^^RITE (6,1031 J“-2,NJ1) 


WRITE(6,104J 
WRITE(6, 105J 
WRI TEI6,106 ) 
WRITE(6, 107) 

102 FORMAT!/, IX, 

103 FORMAn/TlX, 

104 FORMAT!/, IX, 

105 FORMAT!/, IX, 

106 FORMAT!/, IX, 

107 FORMAT!/, IX, 
999 CONTINUE 

RETURN 

END 


CC( J) , 
! D! J J , 
!fc!J) , 
!F( JJ , 

*A* ,ax 

' B* t8X 

*c* ,ax 
^ D* ,ax 

• E* ,3X 
' F* ,8X 


J=2,NJ1) 

J=2,NJl) 

J=2,NJ1) 

J=2,NJ1) 

, IPlOEll. 
, IPlOEll. 
rlPlOEll, 
, LPIOEIU 
, IPIOEU. 
T IP lOEil, 


4/ ! lOX, 
4/! iOX , 
4/( lOX, 
4/!10X, 
4/(1 OX, 
4/( IOX, 


IPlOEll. 4)) 
IPlOEll, 4)) 
IPlOEll. 4)1 
IPlOE 11-4) ) 
lPlOEll-4)) 
IPlOEll. 4) ) 


13. 


f^ROBtJCIBILOT OF THS 
PAOn TS POOR 
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i s FORTRAN IV 360N-F0-479 3-8 THICK OATF 08/20/7‘3 TIME 

j 

j; SUBROUTINE THICK 

I, C***=«= THIS SUBROUTINE CALCULATES DISP. AND MOM. THTKNESSES **:!■.*«** = 

COMMON/CQMl/UfSO) ,T(30) t H t 5 0 J ,RH0 1 5 0 ) , USH t pOJ ,Z(50J ,DZ(50) 

I ^ C0MMaN/CDM2/NJTNJi!NJZ,NJP,NBLf DX 

I i CaMMDN/C0N6/ALPHATBETA,DELTAT GAMMA, RE, VI SCI 50 ) 

1 COMMaN/CaM7/DELlTQEL2THl2,L)ELO 

ji CaMMaN/CaM8/HEIGHT,WIDTH,HSLOPE,WSLQPEr AVWIDH,EFW IOH,OLOA,AHFAC 

; dimension RUH 50) ,RU2I50) 

UINV=1./U{NJ) 

ii RUINV=0.5/RH0tNJ)/UTNJ) 

j" DO 20 J=L,NJ 

RUK J)=RHO{ J1 *U( J)*RUINV 
ii 20 RU2(J) = RUL( J)’^I 1.-UIJ)*UINV) 

fi DELl = {l.-RUi (2)*2.=f=GAMMAJ+0Z{ 2) 

DEL2=RUU2)^2.>«*GAMMA*a .-U<2)*ALPHA»}'UIWV) ^DZiZ) 

?; DO 30 J=3,NJ 

DELl=DELlt( l.-RUl {J)-RUH J“l) )*DZtJ) 

U 30 DEL2=0CL2 + (RU2IJi+RU2[J-Ut*DZ( J) 

r H12-DEL1/DEL2 

j: UPaL=0.995#U(NJ J 

F DO 40 J-2,NJ 

i IF! UtJ J-GT-UPBL) GO TO 4L 

ji 40 CONTINUE 

:i 4i DELO=ZIJ-l)+OZlJ)=i=I (UP3L-UI J-l) )/(UlJ)-U{ J-l) M 

V RETURN 

K END 


I.: 

1 - 


i; 
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MS FORTRAN IV 360N-F0-479 3-8 HPTKKt] DATE 08/20/79 

|; SUBROUTINE HPTRHO 

I; COMMON/COM1/UI50J ,TI50J tH£ 50) ,RHn(50J ,USH(50) ,Z{50) ,0Z( 50) 

P CDMM0N/C0M2/NJtNJ1,NJ2,NJP,NBL,D)£ 

i COMMON/COM4/UCNTR^TCNTR ,rWALL,P,CPINV,KINVTGASK,ZBL 

f PRINV=P*RrNV 

^ DO 223 J=2fNJ 

n J)=H{J)+CPINV 

I 223 RHO( JJ=PRINV/T{ JJ 

RETURN 

END 

5 

I 


■ ! 


I 


I 

I 


I 


I 






riHE 13. 
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3.3 Listing of CIRCULAR Geometry Subroutines. 

Host of the calculations in these subroutines are carried out in terms 
of the radial distance R(J) from the center line. Relations between R(0) 
and Z(J) is given by 

R(J) HEIGHT -Z(J) 


INPUTl Subroutine ; 

This subroutine provides data for the duct shape. 

HEIGHT = Radius of the duct. 

HSLOPE = Axial variation of the radius defined as: 

HEIGHT at x + Ax = HEIGHT at x+HSLOPE-Ax 

ARFAC = 2. for circular ducts 

SIPSIH and ZDZ Subroutines : 

Given U(J), RHO(J), and Z{J)j subroutine SIPSIH calculates SIP{J) and 
SIM(J). It is called at the beginning of each integration step. Subroutine 
ZDZ, called at the end of the integration step, calculates Z(i3) given U(J), 
RHO(J), SIP(J), and SIM(J). Fortran statements in these routines are almost 
direct translations of the formulae given in Appendices A and B of Ref. 7. 

XSAREA Subroutine : 

This subroutine calculates cross-sectional and surface areas of the 
streams. It also calculates cross-sectional area of the duct called OLDA. 
This area is used in the pressure iteration formulae. 
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;]S HJRTaAN IV 360N-F0-479 3-8 


DATE 08/20/79 TIME 


SUBROUTINE INPUTl 

C0MM0N/CQM8/HEIGI-fT,WI DTHfHSLOPE ,WSLQPE , A VWI DH ,EFW lOH^Ol DA,AKFAC 
HEIGHT=W* 

HSLGPE=HB 

ARFACi*2.0 

RETURN 

END 


SIPSIM 


DATE 03/20/79 


34 


JS FORTRAN IV 360N-F0’479 3-8 


TIMF r 


SUBROUTINE SIPSlMtNP) 

COHMQN/COHl/UISOl .nSOJ , H I 5 0) ,RHO 1 5 0 } , USH I 5 0} t 5 OJ r 02 ( 50 | 
CONMON/CUM2/NJtNJlf NJ2,NJPtNBLr DX 
C0MM0N/CQ^13/SIP(50) »S IM ( 50} ,FLQRAT 
C0MP0N/CCM5/AXI50) t AS 15 0) 

COMMON/COM6/A»-PHA,BETATiJELTAtGAMMATRE,VISa5U J 
COMMaN/COM2/DELl,D£L2,H12tOhLO 

C0MM0N/C0M8/HEIGHT* WIDTH,HSLOPE , W SLOPE t AVWI DH »EFW I OH » CL DA ,AKF AC 

DIMENSION RtSOJ ,RR 150), RS 0150} tRUISOJ 

PI = 3. 141593 

PI16TH=PI/16.0 

FLJRAT=0.0 

Rm=HEIGHT 

DO 10 J=2,NJ 

RUI J}=RHO{ JJ«U{J) 

RIJJ=HEIGHT-Z( J) 

RRIJ}=RiJ)*R( J) 

10 RSQ(J3=0.25^IRlJ)4-RCJ-i))t4:2 

C««=SI'S, EXCEPT SIM(2), REDUCED 3Y 4 TO ABSORB 1/4 APPEARING IN SHJ)/4^ 
DO 30 J=2,NJl 

SIPCJ) = PI16TH*{3.=i‘RU( J1 +RU ( J + 1} J * ( RR< J 3 -R SO ( J + 1) J 
30 FLORAT=FLaRAT+ SIP(J) 

SIMt 2) = RU{2)*GAMMA+IRm+Rt2) )*PI*(RUJ-R{2n 
FLQRAT=FLORAT+ S I M { 2 ) .2 5 

DO 40 J=3,NJ1 

SIM( J}=PI16TH*(3.*RU{ J) +RUC J- U ) * t R SO U 1 -RP (J ) ) 

40 FLQRAT=FLORAT+ SIMIJ) 

SIM[NJ) =PIl6TH*C3.’«^RU{NJ)+RUtNJL) ) *RR (N J1 ) *0. 25 
FLURAT=4.0*FL0RAT+4.0=«=S IMINJ ) 

IF( NP.NE. U GC TO 100 
WRITE{6,43) I SIM( J ) , J=3,NJ) 

WRITEt6,441 t SIP(J),J=2,NJU 

43 FORMAT!/,* SIM*,6X , 1 J E 1 L . 3 / ( L OX , 1 OE 1 1 . 3 ) ) 

44 FORMAT!/,* SIP',6X , lOEl 1.3/ ( lOX, lOEl 1 . 3 ) ) 

WR1TE(6,53J FLURAT 

53 FORMAT {/,IOX,‘FLOW RATE= * , Ell. 4) 

100 RETURN 
END 



35 


■JS FORTRAN IV 360N-FJ-479 3-8 2DZ DATE 03/20/79 TIME L3, ^ 

i SUBROUTINE ZDZ 

1 COMHON/COHL/U (50 J tT 1 50J » HI 5 0 ) , RHO ( 50 ) , USH ( 50) » Z 1 5 OJ t OZ I 50 ) - 

r C0MMaN/C0M2/NJtNJl»NJ2f NJP,N8L,0.X 

COMMQN/COM3/SIP150J ,SIM150)tFLORAT ^ 

’ COMHQN/COM6/ALPHAtDETAtDELTAtGAMMArREtVISCf 50) ^ 

: C0MM0N/C0M8/HEIGHT,WIDTH,HSL0PE , WSLOPE , AVWI DH ,EFW lOH » OLDA ,A RF AC ; 

DIMENSION R£50) ,RU(50) i 

• ; RtNJ)=0, 

' 1 PI=3. 141593 i 

i DO 10 J=2,NJ i 

; 10 RUUJ = RHaU)*U£ JJ i 

R£NJl)=t32.s»'lSIPiNJl)+SIMtNJU )/PI/t5.>S'RU(NJl ) + 3.*RU(NJ )) )**0.5 
5 DO 20 N=4,NJ 

J = NJ + 3-N 

* RIJ-U- (l6.=j:5IP( J-1 )/PI/ I 3,=S=RUtJ-n +RU(.I) } + lA.*SIM( J)/PI/{ 

[ L+KU(J-L)H-R(J)*RIJ) )**0.5 

i 20 CONTINUE 

R ( L )= ( S IH 1 2 ) / ( RU t 2 ) *GAMMA*P r ) +R 1 2 ) *R { 2)) ^‘=<'0 . 5 
DO 40 J^2rNJ 

> Zt JJ = RaJ-K(J) 

I 40 DZt J)=Z£ JJ^Zi J-l) 

^ RETURN 

' END 


I 
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; )S FORTRAN IV 360N-FQ-479 3-8 XSARHA OATF 08/ZJ/79 TIME 

i SUaROLTIN£ XSAREAENP) 

; 0UMMDN/C0M1/U£50) iT{50 ) t HI 5 0 } ,RHQ { 50 J ^ OSH ( 50) »Z t 5 3 J rOZ ( 50 ) 

: CQMM0N/C0M2/NJ»NJltNJ2»NJPf NBLtOX 

r C0MM0N/CQM5/AX{ 50) ,AS(50J 

CGMNDN/CDM8/HEIGHTtWIDTH,HSLOPEtWSLOPE, AVwIDHtCFW IJH,OLOA ,APFAC 
' DIMENSION RI50) ,RSQ(50J 

i PI=3. 1-^1593 

PIDX=PI+DX 
; R.{iJ = HEIGHI 

I DO 10 J=2,NJ 

R{J)=HEIGHT-Z{ J ) 

10 KSi.)(J) = 0.25*tR(J)+R( J-1 n=«=*2 

AXI2)=PI« R[2J*(2.*R( l)-R(3J-R{2) ) 

AS( 2) = PIUX’J^R[ U*2. 
ij DO 20 J=3 tNJ1 

; AX(J) = PI^£RSQ( J)-RS' M) 

20 AStJJ=PIDX’S=(RU)+RIJ- 
AX( NJ)=PI=i=RSG{NJ) 

AS(NJ) = PIUX^lR(NJ)+RtNJ D) 

ULUA=PI^HEIGHT«HE IGHT 
Ht IGHT=HE IGHT+HSLGPe*DX 
IFENP.NE.l) GO rn IJO 
WRITE 16,30) I AX £J),J=2»NJ) 

30 FORMAT I/,‘ AX ' , 7X , 1 0 El i - 3/ { 1 OX , L OE 1 1 . 3 ) I 
hKITE(6,A0) (AS t J ) , J = 2,NJ ) 

AOFORMAT(/,' AS ' , 7X , I OF 1 1 , 3/ { 1 OX, iOEl 1 . 3 )) 

100 RETURN 
ENJ 


3.4 Sample Computations— Entry Region Pipe Flow 


DATA: (SI Units) 
MAIN program ; 


C^^DYNAMIC AND THERMaOYNAMI C DATA********^*»***=5=={=*=J=S'=«^**4^*55=x 
UCNTR=15-0 
ZBL=0.01 
TCNTR=400,0 
TWALL=300,0 
P=70927.5 
CP=1003.5 
GASK=i-4 
GASR=287.0 
PRNDL=0.7 
DPOX--5 000.0 

C ❖^GEOMETRIC DUCT DATA**=*=***^^*******its**=^:**4:***i(S:if**+=te*:^!#*5j:*^ 
CALL INPUTl 
X£ND=10.0 

C*<'CONVERGENCE CRI TERI AfGRID SIZEfAND OTHER DA?A’*=********=^==^=* 
DPDTL=0.01 
DXFAC=5.0 
MAXIT=50 
NJ = 21 
NPRNT=10 

C=f'*OATA FOR VSCSTY AND START SUBRDUT I NE *^**’*‘#=«‘ **=«:***« 
VZERO=0. 0000386 
TZ£R0=860.0 
EXPV-.7 
EXPU=7. 

EXPH=7. 

C 


INPUTl Subroutine: 


HEIGHT=0.5 

HSLOPE=0.0 


JOB STATISTICS: 


Machine; IBM 360/44 with DOS operating system. 

Cards Read: 500 

Core Required: 44K 

Compilation Time: 3.19 mns 

CPU Time: 2.11 mns 


VELOUIIY UfTEHPEKATURE T, DENSITY RHOt AND ENTHAPY H DISTRIBUTIONS AS FUNCTION OF Z 
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♦♦OIST. 

along x= 

0.0 

♦♦♦♦NEAR 

WALL RE=l. 

2763E 03 














1 

3. 4106 E- 

-03 

6.8213E-03 

1.0533E-02 

I .60Z0E-02 

2.3519E- 

■02 

3.3 246E- 

-02 

4. 5403F- 

■C2 

6.0179F 

-0? 

7.7751'T 

-02 

9.3291F' 

-02 


1.2L96E 

-01 

1.4892E-0L 

1.7931E-0L 

2. 1329E-01 

2.5098E- 

•01 

2.9254E- 

-01 

3.3610E- 

01 

3.8777E 

-01 

4.9I70E 

-01 

5.0000F 

-01 

U 

1.2863E 

01 

1.4202E 01 

1.50Q0E 01 

1.5000E 

01 

1.5000E 

01 

1.5000E 

01 

1, 5000E 

01 

1.5000E 

01 

1.5900E 

01 

1.5000F 

01 


1. 5OQ0E 

01 

1.5000E 01 

1.5000E 01 

1.5000E 

01 

1.5000E 

01 

1.5000E 

01 

1, 5000E 

01 

1 .inooF 

01 

1.50 OOF 

01 

1 .AonoF 

01 

T 

3.8576E 

02 

3.94686 02 

4.0000E 02 

4,0000E 

02 

4.00COE 

02 

4. OOOOE 

02 

4. OOOOE 

02 

4.0000F 

02 

4.O000F 

02 

&.OOOOF 

02 


4.0000E 

02 

4.0U00E 02 

4.QOOOE 02 

4.U000E 

02 

4. OOOOE 

02 

4. OOOOE 

02 

4. OOOOE 

02 

4. OOOOE 

02 

4.oonoc 

02 

4.0000F 

02 

RHO 

6.4065E- 

-01 

6.,2616E-0l 

6. 1784E-01 

6.1704E- 

-01 

6, 17B4E- 

■01 

6. 178 4E- 

-01 

6. 1784F- 

■01 

6.1784F- 

-01 

6. 1704E- 

-01 

6.l7B4r- 

-01 


6.17fi4E- 

-Q1 

6.17S4E-0i 

6, 17B4E~01 

6 .1784E- 

-01 

6, 17B4F- 

01 

6.17B4F-01 

6, t7B4E- 

•01 

6 .1784E- 

-.01 

6. 17B4E- 

-01 

6.L78■/^F- 

-01. 

H 

3.8711E 

06 

3.9606E 05 

4.0140E 05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.C14QF 

05 

4.014CE 

05 


4.0140E 

05 

4-0140E 05 

4.0140E 05 

4.0140E 

05 

4. 0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 

4.0140E 

05 


+*DIST, ALONG X= 0.2t)00*+**NEAR WALL RE=1.2875E 03 


Z 


3.2813E“03 6.3TeqE-03 1 .0 I95E-02 I . 5A97E-02 2.2629E-02 3.103OE-O2 A.3335E-02 5.7203E-O2 7.3333E-02 9.3117E-02 
1.1525E-QL 1.4031E-01 1.6834E-0L 1.9932E-0L 2.33LOE-01 2.6929E-01 3.0698E-01 3,AA02E-01 3.7511E-OL 5.0000E-OI 


U 


T 


1.3L83E 01 1-A396E 01 
1.580AE 01 1.50CAE 01 


1.5080E 01 1.5582E 01 L.5806E 01 l-'SBOAE 01 1.5804F 01 1.580AE 01 1.9S0AE 01 1.580AE 01 
1.580AE 01 1.5804E 01 1.5804F 01 1.5804C 01 1*5804E 01 1.5004^ 01 l,5*’04=^ 01 1.53n4F Cl 


3.8057E 02 S-aOSlE 02 
3.9990E 02 3.999BE 02 


3.9329E 02 3.9728E 02 3.99a6E 02 4.0000E 02 3.9998E 02 3.9993E 02 3.999BE 02 3.9998E 02 
3.9998E 02 3.9998E 02 3.9998E 02 3.9998E 02 3-9998E 02 3.9998E 02 3.999BE 02 3.9999E 02 


RHO 

H 


6.493DE-C1 6.3604E-01 6.2830C-01 6*2200E-0l £>. 1799E-G1 &.L777E-01 6.1780E-01 6.1779E-01 6. 1780E-01 6.1779E-01 

6.1779E-0L 6.1779E-01 b,1779E-01 6.1779E-01 6. 1779E-01 6.1779E-01 6, 1779E-0L 6.1780E-01 6.1780E-01 6.1778E-01 

3.819QE 05 3.8987E 05 3.9467E 05 3.9867E 05 4.0126E 05 4.0140E 05 4.0L38E 05 4.0138E 05 4.013BE 05 4.013BF 05 

4.O130E 05 4.O130e 05 4.0138E 05 4.0138E 05 4.0133E 05 4.0138E 05 4.0138E 05 4,Q13BE 05 4.0138E 05 4.0139E 05 


**DIST. ALONG X= 1.2000»***NEAB WALL RE=1.3101,E 03 


Z 

u 


T 


3.550.1E-03 7.i089£-03 
1.1633E-01 1.4117E-01 

1. 1853E 01 i.a039E 01 
1.5987E 01 1.5987E 01 


l.Q980E-02 1-.6594E-02 
1.6894E-01 1.9961E-01 

1,3783E 01 1,4462E 01 
1.5987E 01 1.5987E 01 


2.4024E-02 3.3427E-02 
2.3301E-01 2.6873E-01 

1.S025E 01 1.5464E 01 
1.5987E 01 1-5.987E 01 


4.4973E-02 5.B850E-02 
3.0584E-01 3.4215E-01 

1,-577BE 01 1.5959E 01 
l,5987E 01 1.5987E 01 


7.5264E-02 9.4390E-02 
3.723RE-"0l 5.0000E-01 

1.5991E 01 i.59B7E 01 
1,5907E 01 1.5987E 01 


3.7063E 02 3-7790E 02 
3.9990E 02 3.9997E 02 


3.8258E 02 3.8700E ,02 3.909DE 02 3-9425E 02 3.9704E 02 3.9920E 02 4.0007E 02 3.9996E 02 
3.9990E 02 3.9998E 02 3.999BE 02 3.9997E 02 3.9997E 02 3.9997E 02 3.9997E 02 3.9999E 02 


RHO 


6.6670E-01 6.53BSE-01 
6. 1778E-01 6.1779E-01 


6.4588E-01 

6.1779E-01 


6.3849E-01 

6.1779E-01 


6.3212E-QI 
6. L779E-01 


6.2676E-01 

6.1779E-01 


6.2235E-01 
6. 1779E-01 


6.1899E-01 
6, 1780E-01 


8.1764E-01 
6. ITBOE-Ol 


6.1781E-01 

6.1777E-01 


H 


3,7L93E 05 3.7922E 05 3.B392E 05 3.8836E 05 349227E 05 3.9563E 05 3.9843E 05 4.0059E 05 4.0147E 05 4,0136E 05 
4.0138E 05 4,0137E 05 4,0138E 05 4.0138E 05 4.O130E 05 4.0137E 05 4,0137E 05 4,0137E 05 4.0137F 05 4.0139E 05 


♦*□157. ALONG X= 2-2000++++NEAR WALL RE=U3L94E 03 


Z 


3.6739E-03 7.3555E-03 
1.1762E-Q1 1.4236E-01 


1.1354E-02 

1.TOQ3E-01 


1.7139E-02 

2.0058E-01 


2.4775E-02 
2. 3385E-01 


3.4402E-02 

2.6943E-01 


4.6L62E-02 

3.0639E-01 


6,0205F-02 

3.4253F-01 


7.6683E-02 

3.7261E-01 


9.5764F-02 

5.0000E-01 


u 


I. L331E 01 1.248IE 01 
1.6103E 01 1.6100E 01 


1.3214E 01 t,3900£ 01 L.4495E 01 1.4996E 01 1.5407E 01 1.5727E 01 1.5955E 01 1.6085E 01 
L.610LE 01 1.61Q1E 01 1.6101E 01 1.6101E 01 1.6101E 01 1.6101E 01 1.6100E 01 1.610LE 01 


EEPRODUCj^rLlTY OF 'FHE 
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T 

3.6676E 

02 

3. 7367E 

02 

3.78L6E 

02 

i.8245E 

02 

3.8631E 

02 

3.3975E 

02 

3.9277E 

02 

3.9540E 

02 

3.9764B 02 

3.9941E 

02 


4.0005E 

02 

3.9996E 

02 

3.9997E 

02 

3.9997E 

02 

3.9997E 

02 

3.999 7E 

02 

3.9996E 

02 

3.9995E 

02 

3.9996F 02 

3.9999F 

02 

RHQ 

b.7373E- 

01 

ti.6l26E- 

■01 

6-5342E- 

01 

6 .4609E- 

■01 

6.3962E- 

■01 

6.3399E- 

■01 

6.291U-- 

■01 

6.24026- 

-01 

6.?140E-ni 

6.15666- 

01 


6, 1766E- 

01 

6ol781E- 

■OL 

6.17782-01 

6 .1779E- 

•01 

6.1779E- 

•01 

6.1779E- 

■01 

6, 1779F- 

■01 

6.17116 

-01 

6. 17BOF-01 

6.1776E- 

■01 

H 

3.6B04E 

05 

3,74982 

05 

3-7948E 

05 

3.8379E 

05 

3.8767E 

05 

3-9111E 

05 

3.9415E 

05 

3.9679G 

05 

3.9903F 05 

4.0080E 

05 


4.0i45E 

05 

4.0136E 

05 

4.0137E 

05 

4.0137E 

05 

4.0137E 

05 

4.0137E 

05 

4. Cl 36 E 

C5 

4,0135'= 

05 

6.0136C 05 

4.0139E 

05 

=f*=DIilT. 

ALONG X= 3 

i,2000**’*=*NEAR 

WALL RE= 

1. 

3253E 03 













Z 

3.7517E- 

03 

7.5114E- 

■03 

1, 159 IF- 

02 

1 .748 9E- 

02 

2. 5263E- 

02 

3.5049E- 

02 

4.6982E- 

•OZ 

6.11936 

-TJ 

7. 7013F-02 

9.69776- 

■0? 


1.1B82E- 

01 

1.4349E- 

•01 

1,7106E-01 

2.01526- 

■01 

2. 34686- 

■01 

2.7015E- 

01 

3. C6'99E- 

■01 

3. 4302 E- 

-01 

3.7300F-01 

5.0000F- 

01 

U 

i. 1020E 

01 

1.2144E 

01 

1.2866E 

01 

1-3547E 

01 

1,4148E 

01 

1.4669E 

01 

L.SllBE 

01 

1. 548 At 

01 

1. 5783E Cl 

1-6008F 

01 


1.6153E 

01 

1.62Q6E 

01 

1.6201E 

01 

1.6202E 

01 

1. 6-202 E 

01 

1.62026 

01 

1. 6202E- 

01 

1.6202E 

01 

1.6201F 01 

1 .6202E 

01 

T 

3.6431E 

02 

3.7099E 

02 

3.7533E 

02 

3.7951E 

02 

3.8329E 

02 

3.8670E 

02 

3.8975E 

02 

3.9249E 

02 

3.9493E 02 

3.9706F 

02 


3,987gE 

02 

4.00QIE 

02 

3.9996E 

02 

3.9996E 

02 

3.9996E 

02 

3.9996E 

02 

3.9996E 

02 

3.9994E 

02 

3.9995E 02 

3.9999E 

02 

RHO 

6.7824E- 

01 

6.6603E- 

01 

6.5833E- 

■01 

6.5108E- 

•01 

6.446 5E-,0l 

6.3898E- 

■Cl 

6.33976- 

01 

t .2954C 

-01 

6.25(.5F-01 

6.2230E- 

•01 


6. 19&0E- 

01 

6.177 IE- 

■01 

6.1779E- 

01 

6 .1778E- 

•01 

6.1T79E- 

•01 

6.L7T9E- 

•01 

6.170OE- 

■01 

6.17S2E 

-01 

6.1780E-01 

6.L775F- 

■01 

H 

3.6559E 

05 

'S. 7229E 

05 

3.7664E 

05 

3.8084E 

05 

3. 8463 E 

05 

3.8805E 

05 

3.91L2E 

05 

3.93876 

05 

3.9632E 05 

3.9B45E 

05 


4.0018E 

05 

4.0141E 

05 

4,0136E 

05 

4.G136E 

05 

4. 0136E 

05 

4.0136E 

05 

4.0136E 

05 

4.0134E 

05 

4.0135E 05 

A.dlSOE 

05 


**0151. ALCilG X= ^i.200U*#**NEAR V-ALL RE=1.3297E 03 


1 

3.a060E- 

03 

7.6219E- 

■03 

1. 1760E- 

■02 

l,7738E-02 

2.5614E- 

■02 

3.5519E-02 

4.7585E- 

02 

6.1935E- 

■0 2 

7.B690E- 

02 

9.7966E- 

•02 


l. l9asE- 

Ql 

1.445 3E- 

■01 

1.7203E- 

■01 

2.0239E-01 

2. 35*6 E- 

■01 

2,70826-01 

3.0756F- 

■31 

3.43*7r- 

-Oi 

3.7337E- 

01 

5.1000E- 

01 

U 

i.oaosE 

01 

1.1910E 

01 

1.2622E 

01 

1.3298E 01 

I.3899E 

01 

1.4427E 01 

1.4886E 

01 

1.52S1E 

01 

1.5615F 

01 

1.5B39E 

G1 


1.6098E 

01 

1.6239E 

01 

1.6301E 

01 

1.6297E 01 

L.6297E 

01 

1.6297E 01 

1.6297E 

01 

1.6297E 

01 

1.6297E 

01 

1.6297E 

01 

T 

3.62506 

02 

3.690QE 

02 

3,7323E 

02 

3-7731E 02 

3.8102E 

02 

3.8433E G2 

3. 0742E 

02 

3.9018E 

02 

3.9269B 

02 

3.9495E 

02 


3.9696E 

02 

3,9063E 

02 

3.9993E 

02 

3.9996E 02 

3. 9995E 

02 

3.9995E . 02 

3.9995E 

02 

3.9 99 3E 

02 

3.9995E 

02 

3.9998E 

02 

RHQ 

6.8161E-01 

6.6962E- 

•01 

6,6203E- 

01 

6,5487E-0l 

6.4849E-01 

6.4282E-0L 

6.3778F- 

Cl 

6.3327E- 

•01 

6.2922E- 

01 

6.2561F- 

•01 


6.2245E- 

01 

6.1984E- 

■01 

6. 1782E- 

01 

6. 1777E-01 

6. 1779E- 

•01 

6.L779E-01 

6, 17B0E- 

01 

6.1783E- 

■01 

6. IT60E- 

01 

6.1774E- 

01 

H 

3.6377E 

05 

3.7029E 

05 

3.7453E 

05 

3,7863E 05 

3.8235E 

05 

3.8572E 05 

3, B8 77E 

05 

3.9154E 

05 

3.9406E 

05 

3.9633E 

05 


2.9834E 

05 

4e0003£ 

05 

4.Q133E 

05 

4.0136E 05 

4.0I35E 

05 

4.0135E 05 

4.0L35E 

05 

4.0L33E 

05 

4i0136E 

05 

4.0138E 

05 


♦I'DIST. ALONG X= 5.2DOO**!f'*NEAR HALL RE=1.3334E 03 


Z 

3.84916- 

■03 

7.7072E- 

■03 

1.109OE- 

02 

1.7930E- 

■02 

2.5884E- 

■02 

3.5883E- 

■02 

4.8054E- 

■02 

6.2519E- 

-02 

7.93B9E-02 

9.9774E-02 


1.2077E- 

■01 

1.4547E- 

■01 

1.7295E- 

■01 

2.0322E- 

■01 

2- 3619E- 

-01 

2.7146E- 

■01 

3.0809E- 

01 

3.4391F 

-01 

3. 7372 E- 

-01 

5.0000E-D1 

0 

1.0643E 

01 

1.1734E 

01 

1.24376 

01 

1,3108E 

01 

1.3707E 

01 

1.423BE 

01 

1.47 05E 

01 

1.5113E 

01 

1.5468E 

01 

1.5770E 01 


1.6018E 

01 

1,6209E 

01 

1.6337E 

01 

1-6392E 

01 

1.6388E 

01 

1.6389E 

01 

1.63 8 9E 

01 

1.6383E 

01 

1.6388F 

01 

1.63B8E 01 

T 

3.6105E 

02 

3.674QE 

02 

3.7154E 

02 

3.7553E 

02 

3. 7918E 

02 

3.8249E 

02 

3.8550E 

02 

3,88266 

02 

3.9079E 

02 

3.93UE 02 


3.9523E 

02 

3-9712E 

02 

3.9871E 

02 

3.9995E 

02 

3.9995E 

02 

3.999 56 

02 

3.9994E 

02 

3.9992F 

02 

3.9994E 

02 

3.9998E 02 

RHO 

6.B435E” 

01 

6.7252E- 

-01 

6.6503E- 

01 

6.5795E- 

•01 

6.5162E- 

■01 

6.4599E- 

01 

6.4094E- 

01 

6.3639E- 

“01 

6.32276- 

■01 

6.2853E-01 


6.2516E- 

01 

6.2218E- 

•01 

6.1970E-01 

6.1779E- 

•01 

6.1778E-01 

6.1779E- 

01 

6. 1760E- 

01 

6. 17S3E- 

-01 

6, 1780E- 

■01 

6.1773E-01 

H 

3.6231E 

05 

3.6868E 

05 

3.7284E 

05 

3,7685E 

05 

3. 8051E 

05 

3.8383E 

OS 

3. 86 BSE 

C5 

3.8961E 

05 

3.9215E 

05 

3,94496 05 


3.9661E 

OS 

3.9851E 

05 

4.0011E 

05 

4.0135E 

05 

4. 0135E 

05 

4.0L35E 

05 

4. 0134E 

05 

4.0132E 

05 

4.0134E 

05 

4.0138E 05 


*=PDfST. ALONG X= 6.2DaO***#NEAR WALL RE-L.3366E 03 


40 


z 

3.8330E-03 

7,77486-03 

1, 1993E-02 

1,80836- 

-02 

2,60996-02 

3.6173E 

-C2 

4. £4 316- 

■02 

6.29096-02 

7.59586- 

■9? 

9.94406 

-02 


1-2162E-01 

1.4629E-01 

1.7379E-01 

2.O402E' 

-01 

2.36506-01 

2.72U7E-C1 

3. C860E- 

■01 

3.44336-01 

3.74056-01 

5.00006 

-01 

u 

i,0514E 01 

1.1594E 01 

1.2291E 01 

1.29566 

01 

1.3553c Cl 

1.40046 

01 

1.4555E 

01 

1.49726 01 

1.53406 

01 

1 .56606 

01 


1.5933E 01 

1.6153E 01 

1.6329E 01 

1.6441E 

01 

1.6480E G1 

1.64766 

01 

1.64 76E 

Cl 

L.6476F 01 

1. fi4 76E 

01 

1.64766 

01 

T 

3.59826 02 

3.66056 02 

3.7011E 02 

3. 74036 

02 

3.7762E 02 

3.80806 

02 

3.B3 06E 

02 

3.36606 02 

3. B914E 

02 

3.9148E 

02 


3.9366E 02 

3.9S66E 02 

3.9745E 02 

3.98956 

02 

4.0000E 02 

3.9993F 

02 

3.99936 

02 

3.99916 02 

3,99936 

02 

3.99986 

02 

RHO 

6.0667E-O1 

6. 7499E-01 

6. 67586-01 

6 .60586- 

-0 1 

6.54306-01 

6.40706 

-01 

6.43 66 6- 

■01 

^.39106-01 

6.3494E- 

■91 

F..3U3E 

-01 


6.2765E-01 

6. 24486-01 

6.21656-01 

6 .19326 

-01 

6. 17706-01 

6.17806 

-0 1 

6. 1780E- 

01 

6.1 7046-01 

0. 176 JE- 

■01 

6.17726- 

-01 

H 

3.6108E C5 

3.6733E 05 

3.71406 05 

3.75346 

05 

3.78546 05 

3.H222E 

05 

3.85216 

05 

3.3796F OS 

3.90506 

05 

’.9205F 

05 


3.95046 C5 

3.97046 05 

3.98046 05 

4,00346 

05 

4.01406 05 

4.0133E 

05 

4. '91 336 

05 

4.91316 95 

4.01 336 

05 

4.01306 

05 

**DIST, 

ALONG X= 7.-2000+:**’»N6AP 

hALL RE-1.3393E 03 







' 





L 

3.9104E-C3 

7.B297E-03 

1.20776-02 

1.82076- 

-02 

2.62756-02 

3.6410E 

-02 

4. 07396- 

02 

6.33776-02 

0.O43OF- 

■02 

9.9997E- 

-02 


1.22166-01 

1.4699E-01 

1. 74536-01 

2.04776-01 

2. 37586-01 

2.72656- 

-01 

3.09096-01 

3. 44726-01 

3.74306- 

■0,1 

5.00006- 

-01 

U 

1.0410E 01 

1.1479E 01 

1.21706 01 

1.2831E 

01 

1.34266 01 

1.3956E 

01 

1.44306 

01 

1.40526 01 

1.5228F 

01 

1.55616 

01 


1.5852E 01 

1,60996 01 

1.63016 01 

1.6451E 

01 

1.6544E 01 

1.6 5646 

01 

1.65616 

01 

1.65616 01 

1.65616 

01 

1.65616 

01 

T 

3-5876E 02 

3.6407E 02 

3.68876 02 

3.72736 

02 

3.76266 02 

3.79486 

02 

3182426 

02 

3.85146 02 

3.8767E 

02 

3.90036 

oz 


3.9223E 02 

3,94206 02 

3.9617E 02 

3.9707E 

02 

3.992BE 02 

4.00026 

02 

3.99916 

02 

3.9990E 02 

3.9992E 

0? 

3,99986 

02 

RHO 

6.8869E-01 

6. 77156-01 

6 , 6982 E” 0 1 

6 .6288E- 

-0 1 

8. 5666E-01 

6. 51096-01 

o , 46 J 7 E- 

01 

6.4151E-01 

6.37336- 

■01 

6.33486-01 


6.29926-01 

6.266SE-01 

6.23666-01 

6.20996- 

-01 

6. 1B80E-01 

6.17656- 

-01 

6.1782F- 

01 

6.17846-01 

6,17806- 

01 

6.1T71E-01 

H 

3.60016 C5 

3.66156 05 

3.70166 05 

3.74036 

05 

3.77576 05 

3.80806 

05 

3. 83 76 1 

05 

3.3645= 95 

3.0'(O?r 

9 5 

3.91 196 

C' 


3.9360E 05 

3.V566E 05 

3.9756E 05 

3.99266 

05 

4.00606 05 

4.0 14 26 

05 

4. 01 31 t 

05 

4.0130E 05 

4.0! 326 

05 

4.013BE 

05 


**DIST. ALONG X= 8. 2Q00**'!‘*NEAR WALL RE=1.3415E 03 


Z 

3.9331E-C3 

7,87556- 

•03 

1.21476-02 

1.33126-02 

2.64236-02 

3.66106- 

■02 

4. 8999E-02 

6.37046- 

•02 

8.08306- 

02 


1,22716-01 

1.47606- 

■01 

1.7519E- 

■01 

2.0544E- 

-01 

2. 3823 E- 

■01 

2,73236- 

-01 

3.0957E- 

■01 

3.45UE- 

•01 

3. 7469 E-01 

u 

1.0320E 01 

1.1381E 

01 

1.2067E 

01 

1 .2724E 

01 

1.3316E 

01 

I.3847E 

01 

1.43216 

01 

1.47476 

01 

1.5129F 

01 


1,57756 01 

1.6039E 

01 

1.6263E 

01 

1.64426 

01 

1.65706 

01 

1.6640E 

01 

1. 6645E 

01 

1.66436 

01 

1.6643E 

01 

I 

3.5781E 02 

3,6303E 

02 

3,67766 

02 

3,71566 

02 

3.75046 

02 

3.78226 

02 

3.8113E 

02 

3.83836 

02 

3. 8634F 

02 


3,9091E 02 

3,93006 

02 

3. 949 5 e 

02 

3.9674E 

02 

3.98336 

02 

3.9964E 

02 

3.99976 

02 

3.99886 

02 

3.9992E 

02 

RHQ 

6.90516-01 

6,79086- 

■01 

6.7183E- 

01 

6.6495E- 

-01 

6. 56786-01 

6.5324E- 

■01 

6.4825E- 

■01 

6,43696-01 

6.39516- 

■01 


6.32036-01 

6.20686-01 

6.2557E- 

01 

6.22756- 

■01 

6.2027E- 

01 

6.18236- 

•01 

6,1772t- 

01 

6,17866-01 

6.17806- 

■01 

H 

3.5906E 05 

3.65106 

05 

3.6905E 

05 

3.72066 

05 

3.76356 

05 

3.79546 

OS 

3. S247E 

05 

3.05LTF 

05 

3.87696 

05 


3.9228E 05 

3.94376 

05 

3.96336 

05 

3.90136 

05 

3.99726 

05 

4.01046 

05 

4,01376 

05 

4.01286 

05 

4.01326 

05 


**DIST. ALONG X= g,2000«**:S'NEAR WALL RE=1.3438E 03 
£ 


3.9528E-03 7.9147E-03 
1.2318E-0L 1.4813E-01 

1,02446 01 l.i296E 01 
U5703E 01 1.5981E 01 

3.5695E 02 3.6280E 02 
3.8969E 02 3.9180E 02 


i-22D7E-02 1.8401E- 
1.7577E-01 2,0605E- 

1.1979E 01 L.2633E 
1.6221E 01 1.6422E 

3.6675E 02 3.7050E 
3,93796 02 3,9565E 


-02 2,65486-02 3.67B0E' 
-01 Z.3B84E-01 2.73786- 

01 1,32226 01 1,37516 

01 1.6578E 01 1.6683E 

02 3-73946 02 3,77076 
02 3.9734E 02 3.9879E 


■02 4.9221E-02 6.39346-02 B.I172E- 
01 3,iaO3E-0l 3. 45496-01 3.7499E- 

01 1.4227E 01 1.4655E 01 1.5042E 

01 1.6727E 01 1.67236 01 1-6723E 

02 3. 7996E 02 3.82636 02 3.85136 
02 3-9989E 02 3.99886 02 3.9991E 


02 l.008flE-01 
01 5.0000E-01 

01 1.5390E 01 

01 1.6724E 01 

02 3.6747E 02 
02 3.9998E 02 


41 


RHO 6.92i5E-0L 6.6004E-01 6.7365E-Q1 6.6684E-0L 6.6071E-01 6.5521E-01 6.5a24E-0t 6.4570E-ni 6.4152E-0L 6.3763E-01 

6,3400E-01. 6.3059E-01 6.2741E-01 6 .2446E-01 6.2179E-QI 6.I954E-01 6.1782E.-01 6.1734E-01 6,17R1E-3I 6.1769E-0I 

ri i.SaZOE 05 3.6415E 05 3.6B04E 06 3.7iaOE 05 3.7525E 05 3,7539F 05 3.3129F 05 3.3397F 05 3.B647E 05 3.8853F 05 

3.9106E 05 3.9317E 05 3.9516E 05 3.9703E 05 3.9673E 05 4 .0019E 95 4.0129E 05 4.3123F 95 4,0131E 05 4.01JRE 35 

=t=*=DlST. ALONG X-iO. 19994=fr*^NEAR WALL RE=1.345BE C3 

Z 3.9696E-03 r.94B7E-C3 1.2259E-02 1.8478E-02 2,6657E-02 3.0926E-02 4.94UE-02 6.4224E-02 5. 1.467E-32 1.0U3F-C1 

1.23S9E-QI 1.4359E-01 L.762BE-01 2.J659E-01 2.3939E-01 2.7431E-01 3. I049E-01 3.-+5B5F-01 3*7529=-31 5 . 3000F-0 I 

U 1.0178E 01 1.1225E 01 1.1903E 01 1.2553E 01 1.3140E 01 1.366BE 01 1.4143E 01 1.4573E 01 1.4963E 01 1.5117F 01 

1.5637E 01 1.5925E 01 1.617BE 01 1.6396E 01 1.6575E 01 1.6707E Cl 1 .c786l 01 1.5'30ftE 01 U4802E 01 1.6803F. 01 

T 3.5616E 02 3.6201E 02 3.6583E 02 3.6953E 02 3, 7292E 02 3,7602E 02 3.7087E 02 3.3152E O2.3.34Q0C 32 3.3f>i4F 02 

3.3355E 02 3,9067E 02 3-9268E 02 3*94SaE 02 3.9635F 02 3*9795 l 02 3.9925E 02 3.9995E 03 3,9929F 02 3.999<>r. r ,2 

RHO 6.9367E’'01 6.6246E-01 6.7534E-01 4*6O50E-OL 6.6250E-01 6.5704E-01 6.5209E-Q1 A.ATSTF-Ol 6.4339E-01 6.3949E-01 

6.35B4E-01 6.3241E-01 6.2917E-01 6.2614E-01 6.2334E-01 6.2084E-01 6.1d82E-01 6.1773E-01 6,1702H-Ql ft.l768E-0l 

H 3.5741E 05 3.6328E 05 3.6711E 05 3.7082E 05 3. 7423E 05 3.7734E 05 3.8020E 05 3.82fl6E 05 2.9534E 05 3eB769F 05 

3.8991E 05 3.9203E 05 3.9405E 05 3.9596E 05 3.9773E 05 3.9934E 05 4.0064E 05 4-3135E 35 4-0129E 05 4.0138F 05 
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EX 


PREiS 


P-REC9 

WAl^^^EAR 

HEAT FLX^ 

B.L. ThK 

91 SP THK 

KGP 7HK 

SHAPE FAC 

P.O 


7.093E 

04 

0,0 

4.072E-01 . 

3.892E 

03 

1 .018E-02 

1.44BE-03 

1 .141F-03 

1.269E 00 

3.*)U6E- 

-04 

7.093E 

04 

-2.432E-02 

4.255E-01 

3.978E 

03 

1.002E-02 

1.380E-03 

l.lOOE-03 

1.2546 00 

7.812E- 

-04 

7,0«52E 

04 

-4.505E-02 

4.407E-01 

4.04SE 

03 

9.893E-03 

1.331E-03 

1.072E-03 

1.242E 00 

1.56 2E- 

-03 

7.092E 

04 

-6.126E-02 

4.526E-01 

4. lOOE 

03 

9.819E-03 

1.295F-03 

1.051E-03 

1.232F 00 

3. 125E-03 

7.092E 

04 

-7.3B7E-02 

4.6L2E-01 

4. 136E 

03 

9.P03E-C3 

1.P70E-U3 

1.038F-03 

1.223E 00 

6.250E- 

-03 

7.092E 

04 

-8.378E-02 

4 u668E— 01 

4. 156E 

03 

9.B83E-03 

1.255E-03 

1. 0326-33 

1.2166 00 

1.250E- 

-02 

7.092E 

04 

-9.189E-02 

4.692E-01 

4. 157E 

03 

1.132E-02 

1.255E-03 

l.04dE-03 

1.208E 00 

Z.500E- 

-02 

7.Q92E 

04 

-9.820E-02 

4.682E-0L 

4.138E 

03 

1.266P-02 

1.276S-Q3 

1.069E-03 

1. 1946 00 

S.OUOE 

-02 

7.092E 

04 

-1.G36E-01 

4.620E-O1 

4.03 7E 

03 

1 .360L-02 

. 1.333E-03 

1 .138E-03 

1.172E OC 

l.OGOE- 

-Ql 

7.092E 

04 

-1.081E-01 

4.515E-01 

3.978E 

03 

1 .4B3E-C2 

1.454P-U3 

1.280F-03 

1.1356 no 

2.U0UE- 

-01 

7.092E 

04 

-1. 135E-01 

4.324E-01 

3.78BE 

03 

2. 005 E- 02 

1.6936-03 

1.563E-03 

1.083E 00 

3 .OGOE- 

-01 

7.092E 

04 

-1. 180E-0L 

4.ia3E-0l 

3.647E 

03 

2.250E-02 

1.919E-03 

I.B34E-33 

1.0466 00 

4.C00E- 

-01 

7.092E 

04 

-1.216E-01 

. 4.069E-01 

3.534E 

03 

2. 789E-0? 

2.139E-03 

2.095E-C3 

1.0196 00 

S.OjdE- 

-01 

7.092E 

04 

-1.252E-01 

3.974E-01 

■ 3. 44 IE 

03 

3.064E-02 

2.352E-03 

2.355E-03 

9.955E-01 

G.UOOiE- 

-01 

7.092E 

04 

-1.279E-01 

3.894E-01 

3.362E 

03 

3.407E-02 

2, 560E-03 

2.607E-03 

9.819E-01 

7.000E-01 

7-092E 

04' 

-1.306E-01 

3.824E-0L 

3.294E 

03 

3.8661—02 

2.763E-03 

2.054E-03 

9.681E-01 

8.CUOE- 

“01 

7.a92E 

04 

-1.333E-01 

3.762E-01 

3-234E 

03 

4,1456-02 

2.963E-03 

3.099E-03 

9.564E-01 

9.0aOE- 

-01 

7.092E 

04 

-1, 360E-01 

3.7O7E-01 

3, 18 IE 

03 

4.38TE- 02 

3.159E-03 

3.337E-03 

9.4656-01 

l.OOOE 

00 

7.092E 

04 

-1-387E-01 

3.650E-OI 

3. 134E 

03 

4.854E-02 

3.3S2E-03 

3.574E-Q3 

9.379F.-01 

l.iOOE 

GO 

7.092E 

04 

-1.414E-01 

3.613E-01 ' 

3.090E 

03 

5.213E-02 

3.542R-03 

3.808E-33 

9.303F-01 

1-200E 

00 

7.092E 

04 

-l,432E-0l 

• 3.572E-01 

3.051E 

03 

5. 4896-02 

3.730F-03 

4.039E-03 

9.237E-01 

1.300E 

00 

7.092E 

04 

-1,450E-01 

3.534E-01 ■ 

3.015E 

03 

5.735E-02 

3.916H-03 

4. 2 63 E- 33 

9.177F-01 

l.^OOE 

00 

7-092E 

04 

-L.466E-01 

3.499E-01 

2.981E 

03 

6. 045 E- 02 

4. LOlE-03 

4.494F-03 

9. 124E-01 

1- 500E 

00 

7.092E 

04 

-1.486E-01 

3.467E-01 

2.950E 

03 

6.4U8E-02 

4.233E-03 

4.720E-03 

9.075F-01 

1.600E 

00 

7.092E 

04 

-1.505E-0i 

3.437E-01 

2.922E 

03 

6 .a 0 BF -02 

4.463E-03 

4.943E-03 

9.030E-01 

1.700E 

00 

7.092E 

04 

-1.523E-01 

3.409E-01 

2.894E 

03 

7.079E-02 

4. 6426-03 

5 . I P4E-n3 

8. 990F-ni 

l.EUOE 

00 

7.092E 

04 

-1.541E-01 

3.383B-0I 

2. 369E 

03 

r.328E-02 

4. 32JE-03 


g.95.’’£-J] 

i.yaoE 

00 

7.092E 

04 

-i.559E-01 

3.358E-01 

2.845F 

03 

7. 5576- 02 

4.995E-03 

5.602E-03 

8.917E-01 

2. CODE 

00 

7.092E 

04 

-l,577E-0i 

J.335E-01 

2.822E 

03 

7.945E- 02 

5.16 9E-0 3 

5,B18E-33 

3.834E-01 

2 , looe 

00 

7.092E 

04 

-1.595E-01 

3.3I3E-01 

2.S01E 

03 

8.3186-0? 

5. 342F-03 


0^ nS6r_ni 

2.2G0E 

00 

7.092E 

04 

-i,613E-0L 

3.292E-01 

2.7B0E 

03 

6.623E-02 

5.514E-03 

6.248E-03 

a,825E-01 

2.300E 

00 

7.092E 

04 

-1.631E-01 

3.272E-01 

2. 76 IE 

03 

8. 394 E- 02 

5.684E-03 

6.461E-03 

8.799E-01 

Z.'iOOE 

00 

7.092E 

04 

-i.649E-0i 

3.253E-01 

2-742E 

03 

9.143E-02 

5.8S4E-03 

6.672E-03 

8. 7746-01 

2.500E 

00 

7.092E 

04 

-i.66T£-01 . 

3.235E-01 

2.725E 

03 

9.376E-02 

6.023E-03 

6.883E-03 

8.750R-01 

2.600E 

00 

7.092E 

04 

“1.685E-01 

3.218E-01 

2,70BE 

03 

9. 5 92 E- 02 

6.1B9E-03 

7.O92E-03 

8.727E-01 

2.700E 

00 

7.092E 

04 

-1.703E-0L 

3.202E-01 

2.692E 

03 

9. 982 E- 02 

6.356E-03 

7.3016-33 

8, 706E-01 

2.EQQE 

00 

7.092E 

04 

-1.721E-01 

3.187E-01 

2.676E 

03 

i,('33E-ai 

6.S21E-03 • 

7,5086-03 

a,635E-Ql 

2.900E 

00 

7.092E 

04 

-1.739E-01 

3.172E-01 

2.661E 

03 

1.O63E-01 

6.684E-03 

7.714E-03 

8.665F-01 

3.000E 

00 

7.092E 

04 

-1.757E-ai 

3.157E-01 

2.646E 

03 

i .090E-01 

6.348E-03 

7.920E-03 

8.646E-01 

3.100E 

00 

T.092E 

04 

-1.775E-01 

3-143E-01 

2-632E 

03 

1.1166-01 

7, Ollc-03 

8.125E-03 

B.629E-01 

3.200E 

00 

7.092E 

04 

-1.793E-01 

3.130E-01 

2.619E 

03 

1.140E-01 

7.172E-03 

8,3296-03 

3.6116-01 

3.300E 

00 ■ 

7.091E 

04 

-i-QilE~01 

3.118E-01 

2.606E 

03 

i,l62B-0l 

7.334E-03 

8.533E-03 

8.595E-01 

3.400E 

00 

7. 09 IE 

04 

-1.829E-01 

3.105E-01 

2.593E 

03 

1.1B3E-01 

7.494E-03 

8.735E-03 

8. 579E-01 

3.5UOE 

00 

7.091E 

04 

-Ua47E-01 

3.094E-01 

2-381E 

03 

1.2I5E-01 

7.654E-03 

8.937E-03 

8.565E-01 

3.600E 

00 

7.091E 

04 

-l.a65E“01 

3,0aZE-01 

2.570E 

03 

1,2506-01 

7.B13E-03 

9.139E-03 

B.550E-01 

3.700E 

00 

7.09iE 

04 

-1.883E-01 

3.071E-01 

2.558E 

03 

1.281E-01 

7.971E-03 

9.339E-03 

8.535E-01 

3.8U0E 

00 

7.091E 

04 

-L.901E-01 

3.061E-01 

2.547E 

03 

1.309E-01 

a. 129E-03 

9.5 39F-03 

8.522E-01 

3-900E 

00 

7.091E 

04 

-1.919E-01 

3.050E-01 

2.536E 

03 

1.335E-01 

8.287E-03 

9.739E-03 

8.5Q9E-01 

4.000E 

00 

7.09iE 

04 

-1.937E-01 

3.040E'01 

2.526E 

03 

1.360E-01 

8.443E-03 

9.93BE-03 

B.496E-01 

4.100E 

00 

7.091E 

04 

-1.9S5E-01 

3.031E-01 

2.516E 

03 

1.383E-01 

8.599E-03 

1.D14E-0? 

8.494E-01 

4.2UOE 

00 

7.091E 

04 

-1.973E-01 

3.021E-01 

2.506E 

03 

1.405E-01 

6,7556-03 

1,0336-02 

8.472E-01 

4.JOOE 

00 

7.091E 

04 

-1.991E-01 

3.012E-01 

2.496E 

03 

1.425E-01 

a.9l0E-03 

1-053E-02 

e.460F-0l 

4.400E 

00 

7.091E 

04 

-2,Q09E~01 

3.a03E-Gl 

2.486E 

03 

1.445E-01 

9.064E-03 

1.073E-02 

8, 449F-01 

4.500E 

00 

7.091E 

04 

-2.027E-01 

2.995E-01 

2.477E 

03 

1.479E-01 

9.218E-03 

1.092E-02 

8.4386-01 

4.6Q0E 

00 

7.091E 

04 

-2.045E-01 

2.987E-01 

2.4&8E 

03 

1.512E-0I 

9.372E-03 

1.112E-0Z 

0.42BF-Q1 

4.7Q0E 

00 

7.091E 

04 

“2.063E-01 

2.978E-01 

2.459E 

03 

1.541E-01 

9.525E-03 

1.132E-02 

8,417E-01 

4,a00E 

00 

7-091E 

04 

-2.08iE-0I 

2.971E-0I 

2.451E 

03 

1.56BE-01 

9.677E-03 

1.151E-02 

8.4076-01 

4.900E 

00 

7.091E 

04 

-2.099E-01 

2.963E-01 

2.442E 

03 

1.594E-01 

9.S29E-03 

1.170E-02 

8.397E-01 

5.00QE 

00 

7-091E 

04 

-2.U7E-01 

2.955E-01 

2.434E 03 

1.618E-01 

9.981E-03 

1.190E-02 

B,388E-01 

5.100E 

00 

7.091E 

04 

-2, 135E-01 

2.946E-01 

2.426E 

03 

1.641E-DI 

I.O13E-02 

1.209E-02 

8.379E-01 
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5.200E 

00 

7.091E 

04 

-2.153E-0l 

2.941E-01 

2.418E 

03 

1.663E-01 

1.028E-O2 

1.229E-02 

8. 3706-01 

5.300E 

GO 

7-091E 

04 

-2.17lE^0l 

2.934E-01 

2. 41 IE 

03 

1 .634E-01 

1.044E-02 

I .2486-02 

8.3616-01 

5,400E 

00 

7.091E 

04 

-2. ia9E-0l 

2.927E-01 

2.403E 

03 

1.704E-01 

I.Q59F-02 

l.267E-OP 

a.353E-Gl 

5 . 500E 

00 

7-091E 

04 

-2-207E-01 

2,921E-01 

2.396E 

03 

1.723E-01 

1.074E-02 

1.257E-02 

8.344E-01 

5.60CE 

00 

7.09LE 

04 

-Z.225I--01 

2-914E-31 

2.3B8E 

03 

1.790E-01 

1.099E-02 

1.306E-02 

5.336E-01 

5.700E 

00 

7.091E 

04 

-2.243E-01 

2.908E-01 

2.38LE 

03 

1.782E-01 

1.103E-02 

1.325E-02 

8.32RE-01 

5.U00E 

00 

7.091E 

04 

-2.261E-01 

2.902E-01 

2.374E 

03 

1.8L2E-01 

1. 1L3E-32 

1.344E-02 

8.320E-01 

5.900E 

00 

7,091E 

04 

-2.279E-01 

2.894E-01 

2.307E 

03 

1. 6396-01 

L.133E-02 

U363S-02 

8.312E-01 

6.000E 

00 

7.091E 

04 

-2,297E-01 

2.S90E-01 

2. 36 IE 

03 

1.865F-01 

1. 148E-02 

1.3826-02 

8.3056-01 

6. IGOE 

00 

7.091E 

04 

-2.3L5E-01 

2.884E-CI 

2.354E 

03 

1 .989P-C1 

1. 163E-0? 

1.402E-02 

B.297F-01 

6.200E 

00 

7.091E 

04 

-2.333E-01 

2.878E-0l 

2. 348 E 

03 

l,913F.-0l 

1.173F-02 

1. 4216-02 

8.290F-01 

6.3U0E 

00 

7.09LE 

04 

-2.351E-01 

2.873E-01 

2.341E 

03 

1.935E-01 

1.192E-Q2 

1.440E-02 

8,2836-01 

6.400E 

00 

7.091E 

04 

“2.369E-01 

2.867E-01 

2.335E 

03 

1.957E-C1 

1. 2Q7E-02 

U459E-02 

S.275E-01 

6.500E 

00 

7.091E 

04 

-2.38TE**0L 

2.862E-ai 

2.329E 

03 

1.977E-0L 

1.222E-02 

1. 4786-02 

8.269E-01 

6.600E 

00 

7.091E 

04 

-Z.405E-01 

2.a57E'01 

2^ 323E 

03 

1.997E-01 

1.236F-32 

1.497E-02 

. 8. 2626-01 

6-700E 

00 

7.09LE 

04 

-2.423E-0L 

2.352E-01 

2.3L7E 

03 

2.0L6r-01 

1.25LE-Q2 

1.515E-02 

8.255E-0I 

G.aUOE 

DO 

7. 09 IE 

04 

-2.44LE-0.1 

2.-847E-01 

2. 3UE 

03 

2.034E-01 

1.266E-a2 

1.534F-02. 

8.240E-OI 

6.90aE 

00 

T.091E 

04- 

-2.459E-01 

2.842E-01 

2.305E 

03 

2.058F-0L 

1.280E-02 

1.553F-02 

8. 2426-01 

7.GOOE 

00 

7.091E 

04 

-2.477E-01 

2.837E-01 

2-299E 

03 

2.0B9E-Q1 

1.295E-02 

1.572E-02 

8.236E-01 

7.XOOE 

00 

7. 09 IE 

04 

-2.495E-01 

2.833F-0L 

2,294E 

03 

2.117E-Q1 

1.309E-02 

1.59LE-02 

8.229E-01 

7.200E 

00 

7*091E 

04 

-2.514E-01 

2.a28E-0l 

2,2aaE 

03 

2.144E-0L 

1.324E-02 

1. 6106-02 

8.223C-01 

7.300E 

00 

7.091E 

04 

-2.532E-01 

2.623E-01 

2-283E 

03 

2.170E-0L 

1. 338E-02 

1.679E-02 

8. 2176-01 

7.400E 

00 

7.091E 

04 

-2-550E-01. 

2.819E-aL 

2.277E 

03 

2.L94E-0L 

1.3536-02 

1.648E-0? 

8.211E-01 

7.500E 

00 

7.091E 

04 

-2.568E-0L 

2.814E-01 

2.272E 

03 

2.2LTE-0L 

1.367E-02 

1.666E-02 

8.205E-QI 

7.600E 

00 

7.091E 

04 

-2.586E-0L 

2.8L0E-01 

2.267E 

03 

2.240E-0L 

1.332E-02 

1.685E-0? 

8.1996-01 

7*70OE 

00 

7.091E 

04 

-2,604E-0l 

2.B06E-OL 

2*262E 

03 

2.261E-01 

1. 3966-32 

1.704E-02 

a. 1946-01 

7.80QE 

00 

7.091E 

04 

-2.622E-01 

2.802E-01 

2.257E 

03 

2.282E-01 

1.4UE-02 

1.723E-02 

R.iaSE-Ol 

7.900t 

00 

7.091E 

04 

-2.640E-01 

2.798E-01 

2.2928 

03 

2.302E-Q1 

1.425F-02 

1.7426-07 

8.183F-CI 

8.GU0E 

00 

7.091E 

04 

-2.65BE-01 

2.793E-01 

2.247E 

03 

2.321E-01 

1.440E-02 

l,760F-n? 

P. 178F-01 

ti- IGOE 

00 

7.091E 

04 

-2.676E-01 

2. 789E-01 

2.242E 

03 

2.339E-QI 

J.454E-02 

1.779E-02 

3.172E-01 

G.20GE 

00 

7.091E 

04 

-2.694E-01 

2*786E-0t 

2.237E 

03 

2.397E-01 

l.4.6ac-07 

1.798B-02 

a. 1676-01 

6.300E 

00 

7.091E 

04 

-2, 7i2E-01 

2,732E-0L 

2.232E 

03 

2.375E-01 

1.483E-J2 

1.3176-02 

8. 1626-01 

8.40aE 

00 

7.091E 

04 

-2.73aE-01 

2.778E-01 

2.227E 

03 

2.399E-01 

1.4976-02 

1.835E-02 

8.157E-Q1 

G.SUOE 

00 

7.091E 

04 

-2.748E-01 

2.774E-01 

2.223E 

03 

2.428E- 01 

1.511E-02 

1.0 546-02 

8. 1526-01 

8 e 600 E 

00 

7.091E 

04 

-2.766E-0l 

2.771E-01 

2.2i8E 

03 

2.455E-01 

1.526E-02 

1.373B-02 

B.147E-01 

a.700E 

■00 

7.091E 

04 

-2.784E-0i 

2.767E-01 

2.214E 

03 

2.481 E-01 

1.5406-02 

l.89ir-02 

8. 1426-01 

a.sooE 

00 

7.091E 

04 

-z-soaE-Qi 

2.764E-01 

2.209E 

03 

2. 506 E- 01 

1.554E-02 

1,9106-02 

B.137E-01 

a*gooE 

00 

7.091E 

04 

-2.820E-01 

2.760E-01 

2,205E 

03 

2.S30E- 01 

1.56 8E-02 

1.929E-02 

3.132E-01 

9.000E 

00 

7.091E 

04 

-2. 838E-01 

2.757E-01 

2.200E 

03 

2.553E-01 

1. 583E-02 

1.947E-02 

8. 127E-01 

9.100E 

00 

7,a91E 

04 

-2.856E-01 

2.753E-01 

2. 196E 

03 

2.5T4E~01 

1.597E-02 

1,9666-02 

8.122E-01 

9.200E 

00 

7.091E 

04 

-2.B74E-01 

2.7S0E-QI 

2. 192E 

03 

2.596E-01 

1.611E-02 

1.9B5E-02 

R. 1176-01 

9.3UQE 

00 

7.091E 

04 

-2.a92E-01 

2.747E-01 

2*ie7E 

03 

2.616E-0L 

i.625E“02 

2.003E-02 

8.112E-01 

9.400E 

00 

7.091E 

04 

-2.9i0£-0i 

2.743E-01 

2,183E 

03 

2.636E-01 

1.639E-02 

2.022E-02 

0.107E-01 

9.50 0E 

00 

7.091E 

04 

-2.928E-01 

2.740E-01 

2.179E 

03 

2.655E-01 

1.6536-02 

2.041E-02 

a,103E-01 

9.6a0£ 

00 

7-.091E 

04 

-2.946E-0I 

2.737E-01 

2.175E 

03 

2, 673 E-01 

1,6686-02 

2.059F-02 

a. 0986-01 

9*700E 

00 

7.091E 

04 

-2.964E-01 

2,734E-01 

2. 171E 

03 

2-691E~Ql 

l,682E-0? 

2. 0786-02 

R.093E-01 

9,0OOE 

00 

7.091E 

04 

-2.982E-01 

2.731E-01 

2.167E 

03 

2.709E-01 

1.696E-02 

2.097E-02 

a. 0896-01 

9.900E 

00 

7-091E 

04 

-3,OOOE-Ol 

2,728E-0l 

2.163E 

03 

2.726E-01 

1.710E-02 

2.U5E-02 

3.084E-01 

i.QOOE 

01 

7.091E 

04 

-3.010E-01 

2.725E-01 

2.159E 

03 

2.743E-01 

1.724E-02 

2. 134E-02 

8,080E-tU 

i.OiOE 

01 

7.09IE 

04 

-3.036E“0L 

2.722E-0i 

2,155E 

03 

2.771E-01 

1.738E-02 

2.152E-02 

8.0T5E-01 

1.020E 

01 

7.09LE 

04 

-3.054E-01 

2.719E-01 

2.151E 

03 

2.798E-01 

U752E-02 

2.171E-02 

8. 07 IE- 01 
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3.5 Listing of RECTANGULAR Geometry Subroutines. 


For rectangular ducts the computations are carried out along one pair 
of opposing walls. Half the distance between these two walls is called 
HEIGHT. Half the distance between the other pair of opposing walls is 
called WIDTH. The effect of the WIDTH walls is taken into account, as 
in Refs. 5 and 9, by using effective width, EFWIDH, in calculations in- 
volving mass flow rates. Effective width is calculated as, 

EFWIDH = WIDTH - 6^ for Sq < WIDTH 

EFWIDH = WIDTH (l for 6 q> WIDTH 

where 5 q and are the boundary layer thickness and the displacement 
thickness computed along HEIGHT walls. 

It is assumed that the flow is symmetric about the midplanes. 

The computations are carried out from wall to half the height. 

INPUT! Subroutine : 

This subroutine provides data for the duct shape. 

WIDTH =* Half the duct width, 

HEIGHT = Half the duct height, 

WSLOPEjHSLOPE = Slope of the duct walls defined as: 

WIDTH at x + Ax = WIDTH at x + WSLOPE*Ax 
HEIGHT at x+Ax = HEIGHT at x + HSLOPE-Ax 
ARFAC = 1 for rectangular ducts 

SIP5IM and 2DZ Subroutines : 

Given U(0), RHO(J), and Z(J) subroutine SIP5IM calculates SIP(J) 
and SIM(J). It is called at the beginning of each integration step. 


Subroutine ZDZ, called at the end of the integration step, calculates 
Z(d) given U(J), RHO(O), SIP(J), and SIM(O). Width used in these cal- 
culations is the effective width, EFWIDH. Fortran statements in these 
routiens are almost direct translations of the formulae given in Appendices 
A and B of Ref. 7. 

XSAREA Subroutine: 

This subroutine calculates cross-sectional and surface areas of the 
streams. It also calculates the cross-sectional area of the duct called 
OLDA. This area is used in the pressure iteration formulae. 


TV/ 


FORTRAN IV 360N-Fa-479 3-8 


INPUT I 


DATE Od/20/79 


T IME 




SUBROUTINE INPUTl 
CQM.MaN/CaM8/HEIGHTTWlDTH,HSL0PE,WSLDPE, AVWI DH ,EFW IDHi CLOA tA^F AC 
WIDTH=1 
HEIGHT=| 

WSLQPE=1 
HSLOPE^* 

ARFAC=1. 

RETURN 
END 


14 


£ 

i- 


S FORTRAN 
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IV 360N”FO-479 3-8 SIPSIM DATE 08/20/79 TIME 14 


SUBROUTINE SIPSIMINPJ 

CQMNON/COMl/U (50 ) ,T(50) i HI 50) ,RHO(5 0 i ,USH(50) , Z 1 5 0) r DZ { 50 ) 
C0MM0N/C0H2 /NJtNJ1tNJ2i NJP» NQLtOX 
COMMON/COM3/SIPI50) ,SIM(50) ,FLORAT 
CQMM0N/CCJM5/AX(50)tAS (50) 

C0MHDN/CDM6/flLPHA,BETA,DELTA,GAMMA,RETVISC( 50) 
C0MMGN/CaM7/0ELl,DEL2,HI2,0EL0 

COMMON/COMa/HEIGHr, MI DTH tHSLDPE t WSLDPE, AVM lDH,EFW U)H, GLOA,ARFAC 
dimension RU(50) 

If IOLDA.L£.O.Q)EFWIOH=N lOTH-DELl 
WIDTT^EFWIOH^O. 03125 
FLGRAT=0.0 

C=!==i'SI ‘StEXCEPT SrM£2)f REDUCED BY 4 TO ASSDRV 1/4 APPEARING IN Sl(J)/4«# 
DO 10 J=2,NJ 
10 RU( J I^RHOtJ I’S'UI J) 

DO 30 J-2rNJl 

SIP( J) = I3.>!'RU( JI+RUI Jfl ) )«0Z( J+imJIDTT 
30 FLORAT=f LCRAT+SIP ( J) 

SIM(2) = RU(2)*GAMMA*DZ(2 J*WIDTT=^32. 

FLURAT=FLORAT+S IM (2 )=J=.2 5 
DO 40 J=3,NJ 

SIM(J) = {3,*RU{J1 ^-RU^ J-1) )*DZ( J) *WIDTT 
40 FLURAT=FLORAT+S I.M(J ) 

FLURAT=fLCRAT={=4,0 
IF (NP.NE.l) GQ TO 100 
WRITE(6t22) ( SIM(J) , J=2tNJ) 

WRITE(6,24) ( SIP(J), J=2tNJl) 

WRITE(6,53) FLCRAT 

22 FORMATf/,' , lOEil . 3 / (1 OX 1 1 OEl 1 . 3 )) 

24 FORMAT!/,^ SIP'»6X t 10 El 1 . 3/ { 1 OX , I OE 1 1 . 3 ) ] 

53 FORMAT £lX,9Xt*FLCw RATE= • , Eil.Aj 
100 RETURN 
END 
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FORTRAN IV 360N-FJ-479 3-8 ZOZ DATF 08/20/79 TIMF 14 

SUBROUTINE ZD2 

CONNCN/C0Ml/U(50) tT(50) ,H{50} ,RHO ( 50 ] ,USH ( 50) , Z I 5 0) , DZ < 50 ) 
CQMM0N/C0N2 /NJtNJ1,NJ2,NJPtNBL,DX 
COMMQN/COH3/SIP(50) ,S IM { 50 } » FLORAT 
CGMMDN/C0M6/ALPHA,BFTA,DELTAtGAMMAtRE,VISC(50 J 

CDMMDN/CQM8/HElGHTr WIDTH, HS LOPE ,WSLOPE , AVWI OH ,EFW lOHiOLDA ,AKF AC 
DIMENSION RU150) 

WINV= 8.0/EFWlDH 
DO LO J=2,NJ 
10 RU( J)=RHDIJ )«U{ J) 

ZJ2)=.125*WINV*SrM( 2)/RUI2)/GAMHA 
D2{2)=Z(2) 

DO 20 J=3,NJ 

DZ (J ) = ( S I P i J- 1 J +S IM I J ) J / ( RU ( J ) +RU (J -1 n -W IN V 
20 ZI J) = Z( J-1) ^-DZ^J) 

RETURN 

END 



f;IS FORTRAN IV 360N-F0-A79 3-8 XSAREA 
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DATE 03/20/79 TIMC 14 


SUBROUTINt XSAREAINPJ 

C0MH0N/C0M1/U150) »T ( 50J , HI SO J ,RHO ( 50 ) , USH t 5 0) tZ ( 5 0) fOZI 50) 

COMWQN/COM2/NJiNJIrNJ2rNJPtNBLTOX 

CDMMQN/CGM5/AXI50) f ASI50J 

C0MH0N/C0M7/D£Ll,DEL2,H12iUEL0 

C0MK0N/C0H8/HEIGHTt WIDTH tHSLOPE ,WSLOPE, AVWIOH ,EFW lOHiOLOA tARF AC 
C^=^CALCULATION OF NEW OUCT CROSS-SECT tQN4'=4=*=5'RRR RRRRRRRRRRRRRR RR RRRR PRRRR 
HE! GHT=HEIGHT+HSLCPE=<'OX 
WIDTH=WIDTH+WSLDPE=f=DX 
AVWlDH=WIDTH + WSL0PE=:'DX=i'0,5 
WIUC0R=DEL1 

IF! OELO.GT.WIOTH) WIOCOR =W10TH*0EL 1/DElO 

EFWIDH=HIOTH-WIDCOR 

ULOA =HEIGHT*EFW1DH 

WHF-0 .S^^AVW IDH 


AXI2)=WHF*IZI3J+ZI2)) 

DO 10 J=3tHvJ1 

10 AXI J)=WHF#( Z( J+U-Zt J-1 ) ) 

AX( NJ)=WHF« {ZINJ)-ZINJ1 J ) 

00 20 J=liNJ 
20 ASl J)=DX*AVWIDH 

IF! NP.NE, UGCl TO 35 
WRITEI6,30) I AS! J) r J = 2,NJ) 

30 FUR^^ATC/,' AS* ,7X,10E11 .3/1 IOXt lOEll.3) 1 
WRITE (6,32) (AX{ J) ,J=2,NJ) 

32 FORMAT I/r* AX * , 7 X , 10 El 1 . 3/ (1 CX , IGE 1 1 . 3 ) ) 
35 RETURN 


END 



i 


3.6 Sample Computations— Rectangular Diffuser Flow 


DATA: (SI Units) 


MAIN Program : 


C^^DYNAHlC AND THERMODYNAMIC 
UCNTR=200.0 
ZBL=0,015 
TCNTR=300.0 
TWALL=300*0 
P=71000.0 
CP=1003.5 
GASK=1,4 
GASR=28T.O 
PRNDL=0,7 
DPDX=210000.0 
C**GEOHETRIC DUCT 
CALL INPUTl 
X£ND=0.5 

C**CQNVERGENCE CRITERIAf GRI D SIZE»ANO OTHER DAT 
0PDTL==0-01 
DXFAC=5.0 
MAXIT=50 
NJ=21 
NPRNT=5 

C*«DATA FDR VSCSTY AND START S UBR OUT I NES* 
VZER0=0,0000386 
TZERO=860.0 
EXPV=-7 
EXPU=7, 

EXPH=7, 

C 


INPUTl Subroutine : 

WIJTH=0.025 

HEIGHT=0,05 

WSLOPE=0.03 

HSL0PE=0,03 

JOB STATISTICS : 

Machine: IBM 360/44 with DOS operating system. 

Cards Read: 491 

Core Required: 42 K 

Compilation Time; 3.05 mns 

CPU Time: 1.97 mns 
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VELUCITY U, TEMPERATURE T, DENSITY RHOt AND ENThAPY H DI STRIBUTtONS AS FUNCTIC'i OF I, 


* 


j **DIST. 

ALONG X= 0,0 

**««NEAR 

WALL RE=1. 

200GE 03 







! 1 
5 

2.A227E-0A 

l.lOOHE-02 

4.8453E~04 

1.3637E-02 

7, 7459E-04 
1.6648E-02 

1.2226E-03 
2 .Q057E-02 

1.8568E-03 

Z,3aa4E-02 

2.7 J40E-C3 
2.B151E-02 

3. 7399C-0^ 
3.2376E-02 

5.1370I;-T2 

3,8079E-0? 

77546-03 

4.3781E-0? 

8.7222F-03 

=:.O000E-02 

t 

• u 

1,I093E 02 
l.9134e 02 

1.2248E 02 
1.973a£ 02 

1.3097E 02 
2.0000E 02 

1.3979E 02 
2.0000E 02 

1.4839E 02 
2.0000E 02 

1.5658E 02 
2.OOU0E 02 

1.6432E 0? 
2.0QOOE 02 

1.7162E 02 
2.00J0E 0? 

I.T553F 3? 
2.0000E '?2 

1.35D9F 0? 
2..1000F 02 

' T 

I 

3.0000E 02 
3,oooqE 02 

3.000CE 02 
3.JOOOE 02 

3.0000E 02 
3.0000E 02 

3.0COOE 02 
3.0000E 02 

3,OOOOE 02 
3.0000E 02 

3.0000E 02 
3.0000E 02 

3.C000E 02 
3.0000E 02 

3.0000F 02 
3.7 001)6 3 2 

3.0000E 02 
3.0300E 02 

3.0000E 
3,nOOOE 02 

i RHO 

8,2462E-01 

8.2462E-01 

8.2462E-01 

8.2462E~01 

8.2462E-01 
8.2462 E-01 

d.2462E-01 

8.2462E-01 

8.2462E-01 

8.2462E-01 

8.2462E-01 

8.2462E-0L 

3.2462E-0L 

8.2462E-01 

8.2462E-JI 

ri.2462'=-0l 

8.2462E-01 
3. 2462F-01 

8.2462E-01 

3.2462F-91 

H ■ 

3,Q105E 05 
3.0105E 05 

3.0105E 05 
3.0105E 05 

3.0105E 05 
3.D105E 05 

3.01D5E 05 
3.0105E 05 

3.0105E 05 
3.0105E 05 

3,0105E 05 
3.0i05E 05 

3. OUSE 05 
3. 01 05 E 05 

3.0105F 05 
3.0105E 05 

3.01J'5E. 05 
3.0105E 05 

3.0105F 05 
3.0105E 05 

♦vOIST- 

ALONG X= 0.0006****‘NEAR 

HALL RE=i. 

iOaOE 03 







l 

2.4346E-04 

l,10llE-02 

4,3646E-04 

1.3647E-02 

T..76S9E-04 

1.6660E-02 

1.2254E-03 

2.0070E-02 

1.86C3E-03 
2. 3898E-02 

2.70836-03 
3.816 5E-02 

3.7951E-03 

3.2991F-02 

5.1452E-03 

3.3096E-02 

6.7325E-03 

4.379«E-02 

8.7303E-03 

5.00I9F-02 

U 

1.L034C 02 
i.9109E 02 

1,2219E 02 
1.9696E 02 

I.306SE 02 
L.9974E 02 

1.3949E 02 
L.9975E 02 

1.48C9R 02 
1.9975E 02 

1.5628E 02 
1.9975E 02 

1.6402E 02 
1.9975E 02 

1.7134E 02 
1.9975E 02 

1.7826E 02 
U9975E 02 

1,84B2E C2 
1.9975E 02 

1 

! T 

J.0020E 02 
3.0005E 02 

3,0aC5E 02 
3.0005E 02 

3.0J05E 02 
3.0005E 02 

j.OJJEE 02 
3.0004E 02 

3.0e05E 02 
3.0005E 02 

3.000SL- 02 
3.0004E 02 

3.C005L 02 
3.0004E 02 

3.0005E 32 
3.0004E 02 

3.0005F )2 

3. 0004>= 02 

3.0005E C2 
3.f)005F 02 

RHO 

8.2450E-01 

8.2492E-0L 

8.2491E-01 

B.2492E-0i 

0.2492E-C1 

8.2492E-01 

8.2492E-9I 
3 , 2493E— 0 1 

B.2492E-0I 

a,2493E-0l 

a.2492E-0l 

8.2493E-01 

e.2492r-oi 
8. 24936-0 i 

8.24926-0 1 
8.2493F-01 

“. 24 02E-01 
8.2493E-QI 

H..»4926- ;1 
8.2492E-01 

i H 

1 

) 

3.0125E 05 
3.0110E 05 

3.0U0E 05 
3,OllOE 05 

3.0110E 05 
3.0110E 05 

3.Q11QE 05 
3,0109E 05 

3.01100 05 
3.0110E 05 

3.0HOE 05 
3.0109E 05 

3.0110E 05 
3.0109E 05 

3.0110E 05 
3.0109E 05 

3.0110E 05 
3.0109F 05 

3.0110E 05 
3.01106 05 

*wDisr. 

ALONG X= 0.0200**>i'*NEAR 

HALL RE=1. 

1707E 03 







z 

2.6072E-04 

I.1277E-02 

5.1B72E-04 

1.3953E-02 

8.2382E-04 

1.7000E-02 

1.2912E-03 

2.0448F--02 

1.9485E-03 

2.4301E-02 

2.8220E-03 

2.85986-02 

3.9367E-03 
3.33 56E-02 

5.3165^-33 

3.8597F-02 

6.9844E-03 

4.4339E-0? 

B.9632F-03 

5.0601E-02 

i 

1.0061E 02 
l.a289E 02 

1.1271E 02 
I.BSOaE 02 

1.2151E 02 
1.9174E 02 

1.3054E 02 
1.9268E 02 

1.3931E 02 
U9255E 02 

1.476 5£ 02 
1.9257E 02 

1.5555E 02 
1-9257E 02 

1.6303E 02 
1.9257E 02 

1.7013E 02 
I.9257E 02 

1.7680E C2 
1.9257E 02 

r 

3.02B2E 02 
3.0147E 02 

3.0235E 02 
3.(U46E 02 

3.0204E 02 
3.0145E 02 

3.0179E 02 
3.0144E 02 

3.0163E 02 
3.0144E 02 

3.0155E 02 
3.0144E 02 

3.0152B 02 
3.01446 02 

3.0150E 02 
3.0144E 02 

3.0149E 02 
3.0L44E 02 

3.0148E 02 
3.0146E 0? 

; RHO 

8.2966E-01 

0.3338E-O1 

8.3097E-01 

8.3341E-0L 

8.3182^01 

8.3343E-01 

8.3251E-01 
8 ,3346E-01 

a. 32 94 E- 01 
3. 3345E-01 

P.3315E-01 

8.3346E-01 

8.3325E-Q1 

8.3345E-01 

3. 3330E-01 
8.3345E-01 

3. 3333 E- 01 
a.3345E-0l 

B.3336E-01 

8.3342E-01 

H 

3.Q388E 05 
3.0253E 05 

3.0341E 05 
3.0252E 05 

3.0309E 05 
3.0251E 05 

3.0284E 05 
3.0250E 05 

3, 0269E 05 
3,02506 05 

3.0261E 05 
3.0250E 05 

3.0257E 05 
3.0250E 05 

3.0256E 05 
3.0250E 05 

3.0254E 05 
3.0250E 05 

3.0253E 05 
3.0251E 05 

#*DIST. 

ALONG X= 0.0700****fJEAR 

HALL RE=l.i247E 03 







2 

2,8223E~04 

L.1893E-02 

5.6052E-04 

1,4676£-02 

6.8928E-‘04 

1.7839E-02 

1.3916E-03 

2,l375E-02 

2.0953E-03 

2.5307E-02 

3.0258E-03 

2.9683E-02 

4.2071E-03 

3.4532E-02 

5.6617E-03 

3.9871E-02 

7.4121E-03 

4.5721F-02 

9.4810F-03 

5.21026-02 

U 

f 

8.8177E 01 
1,6616E 02 

9.8914E 01 
1.7134E 02 

1.0689E 02 
1,7548E 02 

1.1524E 02 
1.7819E 02 

1.2352E 02 
1.7BB2E 02 

1.3156E 02 
1.7873E 02 

1.393LE 02 
i.7074£ 02 

L.4670E 02 
U7874E 02 

1.5369E 02 
1.7874E 02 

1.6022E 02 
1.7874E 02 


( 
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T 

S.OSi-tZ 02 

3.0519E 

02 

3,0503E 02 

3.0485E 

02 

3.0467E 

02 

3.0450E -02 

3.0437E 02 

3.0426E 

02 

3.0415r 

02 

3.0412F 02 


3.0408£ 02 

3.0405E 

02 

3.0402E 02 

3.0400E 

02 

3.0400E 

02 

3,0400E 02 

3.0400E 02 

3.0400E 

02 

3.0400E 

0? 

3. 0402 F -02 

RHO 

8.45Q3e-01 

8.4545E- 

-01 

8.45B9E-01 

8.4640E- 

-01 

3.4690E- 

■01 

8.4736E-01 

a.4774E-01 

a .4H03E^ 

-01 

0.AR25F- 

■01 

8.4842E-01 


8.4854E-01 

8.4863E- 

-01 

3.4870E-01 

8.4875E- 

-QL 

a. 48 78E- 

■01 

6.4877E-0i 

8.4377E-01 

3.4877E' 

-01 

8. A577f- 

;01 

8.48728-01 

H 

3.0641E 05 

3.U626E 

05 

3.0610E 05 

3.0592E 

05 

3.0573E 

05 

3.0557E 05 

3-0543E 05 

3.0533E 

05 

3.0525E 

05 

3.0519E 05 


3.0515E 05 

3.0511E 

05 

3.0509E 05 

3.0507E 

05 

3.0506E 

05 

3.0506E 05 

3, 0506E 05 

3.0506E 

05 

3.0506F 

05 

3.0508F 05 

't=»Disr. 

ALONG X= 0.1200«‘**+NEAR 

WALL RE=l,i 

0854E 03 











z 

3.0085E-C4 

5.9649E- 

-04 

9.^549E-04 

i.478ZE- 

-03 

2.2Z35E- 

■03 

3.2072E-03 

4. 4530E-03 

5.9833E- 

-93 

7. 8197E- 

■03 

9.9B40C-03 


1.2499E-02 

1.539IE- 

■02 

l.a867E-02 

2*2309E- 

-02 

2.6333E- 

02 

3.078 l)E-02 

3.5720E-02 

4.1 153E- 

-02 

4. no5t- 

02 

5.3597E-02 

u 

7.8898E 01 

8.8637E 

01 

9«5890E 01 

1.03S5E 

02 

1. U21E 

02 

1.1872E 02 

1.260ZE 0? 

1.3306E 

02 

1.39788 

02 

1.4613F: 02 


1.5202E 02 

i,573IE 

02 

1.6179E 02 

1.65X7E 

02 

1.6716E 

02 

1.6733E- 02 

1.67306 02 

1.6731E 

0? 

1.6730E 

02 

l,6730E 02 

T 

3.0684E 02 

3.0692E 

02 

3.0692E 02 

3.Q6B6E 

02 

3- 06 76 E 

02 

3.0665E 02 

3. 0653E 02 

3. 064 IE 

02 

3.0631E 

02 

3.0621E 02 


3.0614E 02 

3,0607E 

02 

3.0G02E 02 

3 .0599E 

02 

3.0597E 

02 

3.0596E 02 

3.D596E 02 

3.0596E 

02 

3.0596E 

02 

3.0599E 02 

RHO 

B,5793E-01 

8*5770E- 

-01 

8.5771E-01 

3.57B7E- 

-01 

Q.5814E- 

01 

8.5846E-01 

8.5880E-01 

8.5912E- 

-01 

8.5942E- 

01 

8.5968E-01 


8.5990E-Q1 

0.6OO7E- 

-01 

8.6021E-01 

8.6031E- 

-01 

0.-6O37E- 

01 

8.6039E-01 

8.6038E-0i 

a.603BE- 

-01 

8.6035E- 

■01 

8.6031E-01 

H 

3.0791E 05 

3.0799E 

05 

3,0799E 05 

3.0793E 

05 

3.0784E 

05 

3.0772E 05 

3.07 60E 05 

3.0748E 

05 

3.0733E 

05 

3.07Z9E 05 


3.0721E 05 

3.0714E 

05 

3.0710E 05 

3.0706E 

05 

3.0704E 

05 

3.0703E 05 

3.0704E 05 

3.0704E 

05 

3.0704E 

05 

3.Q706E 05 


*#Dlsr. ALONG A= 0. 1700***«NEAR WALL HE=l,05L3E 03 


2 

3.1814E- 

■04 

6.2990E- 

■04 

9.9785E- 

■04 

1 .5592E-03 

2. 3439E- 

-03 

3-37B5E- 

-03 

4.68 70E- 

-03 

6,291BE- 

-0 3 

8.2139E-03 

1.0474E 

-02 


i.3095E-02 

1.6099E- 

-02 

1.9490E- 

■02 

2.3245E’02 

2.7373E- 

•02 

3.1915E^ 

-02 

3.692BE-02 

4.Z449E- 

-02 

4.B495E- 

■02 

5.5097E- 

-32 

U 

7.1570E 

01 

8.0455E 

01 

8.7108E 

01 

J.4136E ul 

1.0125E 

02 

1.0824E 

02 

1.150BE 

02 

1.2173E 

02 

U2H14F 

02 

1.3427F 

02 


1.4005E 

02 

1.4537E 

02 

1.5005E 

02 

1.53B7E 02 

1.5656E 

02 

1.5781E 

02 

1.5773E 

02 

1.5774E 

02 

1.5774E 

02 

1.5774E 

02 

T 

3.0787E 

02 

3.0813E 

02 

3.0824E 

02 

3-0829E 02 

3.0828E 

02 

3.0823E 

02 

3.0815E 

02 

3.0806E 

02 

3.0796E 

02 

3.0786E 

02 


3.0777E 

02 

3.0769E 

02 

3.0762E 

02 

3.0756E 02 

3.0753E 

02 

3-075 IE 

02 

3.0751E 

02 

3.0751E 

02 

3.0751E 

02 

3.Q754E 

OH 

RHO 

8.6830E- 

■01 

a.6758E- 

-01 

8.6727E- 

01 

3.6713E-01 

8.6715E- 

01 

8.6730E- 

-01 

6.6751E- 

01 

8.6777E- 

-01 

8.6805E- 

■01 

0.6832E-01 


a.6858E- 

■01 

8.6881E- 

01 

8.6901E- 

01 

8.6917E-01 

8.6927E- 

■01 

8.693 3E- 

-01 

8.6933E- 

01 

8.693ZE- 

-01 

a,6932E- 

•01 

8.69Z4E- 

-01 

H 

3.0895E 

05 

3.092LE 

05 

3.0932E 

05 

3.0937E 05 

3.0936E 

05 

3.093 IE 

05 

3.09Z3E 

os- 

3.0914E 

05 

3.0904E 

05 

3.0R94E 

05 


a.oaasE 

OS' 

3-0877E 

05 

3-.0870E 

05 

3.0864E 05 

3. 0660E 

05 

3.0858E 

05 

3.0858E 

os 

3.oa50E 

05 

3.0858E 

05 

3,0861E 

05 


#*DiSr. ALONG X= 0.2200++**NEAR WALL RE=l-02UE 03 


Z 

3.3490E- 

•04 

6-6236E- 

■04 

1.D486E-03 

1.6377E- 

-03 

2.4606E- 

•03 

3.S448E- 

■03 

4.9146E- 

■03 

6.59246- 

-0 3 

8.5990E- 

■03 

1.0955E- 

■02 


1.3680E-02 

1.6797E- 

•02 

2.0305E-02 

2.4178E- 

-02 

2. 8417E- 

•02 

3.3057E- 

■02 

3. 8153 E- 

■02 

4.3 756E- 

-02 

4.90976- 

■02 

5.6596E- 

■02 

u 

6.5502E 

01 

7-36856 

01 

7-9825E 

01 

8.6352E 

01 

9.2941E 

Ql 

9.9472E 

01 

1.0589E 

02 

1.121TE 

02 

l- 1828E 

02 

1.2418E 

02 


1.2981E 

02 

1.3508E 

02 

1.3985E 

02 

1 .4393E 

02 

1.47UE 

02 

1.491 BE 

02 

1.4966E 

02 

1.4959E 

02 

1.4960E 

02 

1.4959E 

02 

T 

3.0863E 

02 

3.0902E 

02 

3.0922E 

02 

3.0936E 

02 

3.0942E 

02 

3.094 3E 

02 

3.Q940F 

02 

3.0934E 

02 

3.0926F 

02 

3.0910F 

02 


3.0909E 

02 

3.0^006 

02 

3.0892E 

02 

3.O085E 

02 

3.O80OE 

02 

3.0876E 

02 

3.0875E 

02 

3.0875E 

02 

3.08756 

02 

3.0879E 

02 

RHO 

8.76T0E- 

•01 

8.7559E- 

■01 

8.7502E- 

■01 

8.7464E- 

-01 

a. 7446E- 

■01 

8.7443E- 

■01 

8,7452 E- 

■01 

8.74696- 

-01 

8.7490E- 

•01 

8.7514E- 

•01 


8.7539E- 

01 

B.7564E- 

■01 

8,75876-01 

8.7607E- 

-01 

8.7623E- 

-01 

8.7633E- 

■01 

8.7637E- 

01 

B.7635E- 

-01 

8. 7635 E- 

■01 

8.7626E- 

01 

H 

3*09716 

05 

3*1010E 

05 

3.1030E 

05 

3.1044E 

05 

3. 1050E 

05 

3.1051E 

05 

3.1048E 

05 

3.1042E 

05 

3. 1Q35E 

05 

3.1026E 

05 


3.1017E 

05 

3.1008E 

05 

3.1000E 

05 

3.0993E 

05 

3.09886 

05 

3.C984E 

05 

3. 09 83 E 

05 

3.09B3E 

05 

3.0983E 

05 

3.0987E 

05 


^^DIST. ALONG X= 0,2 700^*«'#NEAR VtALL RE=9.fl469E 02 


53 


z 

3.6020E-04 

7.10506-04 

1.1225E-03 

1.74376-03 

2.62C4E-03 

3.7648B-03 

5.2060E- 

■03 

6.«>ft5 3E 

-03 

9.0643E- 

-03 

1.1«2LE-C2 


1.4354E-02 

1.7585E-02 

2.12096-02 

2.51966-02 

2,95436-02 

3.4280E- 

-02 

3. 94526-02 

4.5120E 

-02 

5.1330F- 

-02 

5.8105F-02 

u 

5.83226 01 

6.5B47E 01 

7.15756 01 

7.77206 

01 

8.39656 31 

9,0184E 

01 

9.63256 

01 

1.0236E 

02 

1.0926E 

02 

1.1400P 3? 


1.1952E 02 

1.24776 02 

1.2962E 02 

1.3388E 

02 

1.37416 02 

1.4000F 

02 

1.4138E 

02 

1.41376 

32 

1.4137E 

02 

1.41366 02 

T 

3.09396 02 

3.09896 02 

3.1016E 02 

3.10366 

02 

3-1048E 02 

3.1054E 

02 

3-10556 

02 

3-1053E 

02 

3. 1048E 

02 

3.1041F 02 


-3.10336 02 

3.10256 02 

3.L016E 02 

3.10096 

02 

3,10026 02 

3.Q997E 

02 

3.09946 

02 

3.0994E 

02 

3.0994E 

02 

3.0995C 0? 

KHQ 

8.8454E-01 

8.33116-01 

8. 82346-01 

8.B177E-01 

8.0L42E-O1 

8,81256- 

■01 

8. 8122E- 

•01 

9. 3 1296 

-01 

B.8143E- 

-01 

3.81626-01 


a. 81846-01 

8.32086-01 

8.82326-01 

a. 82546-01 

8.82736-01 

8. 82976-01 

B.8296E- 

'01 

9.8296E 

-01 

0.8296F- 

■01 

8.8284E-C1 

H 

3.1047E 05 

3.10976 05 

3.11246 05 

3.11446 

05 

3-11576 C5 

3.U636 

05 

3.11646 

05 

3.11616 

05 

3. 11 56F 

C5 

3.1150E 05 


3.11426 05 

3.U33E 05 

3. IIZ5E 05 

3,11176 

05 

3. HUE 95 

3. HOSE 

05 

3. 11026 

05 

3.11026 

OS 

J.1102E 

05 

3.1107E 05 

**OIST. 

ALONG 0.320e*++*N£AR 

WALL RE=9.i 

54 09E 02' 












Z 

3.8037E-04 

7.49066-04 

1.18246-03 

1 .8406E-03 

2.75596-03 

3.95596- 

-03 

5.4648 6- 

'03 

7.30366- 

-03 

9.4917E- 

■03 

l,2047F-02 


1.4989E-Q2 

1.83346-02 

2.20766-02 

2 .61816- 

-02 

3. 06406-02 

3.54826-02 

4. 0742 6- 

■02 

4.6479E- 

-02 

5.27546- 

■02 

5.9601E-9? 

U 

5.31946 01 

6.013LE 01 

6.54066 01 

7.10926 

01 

7.69C8E 01 

a. 2 7426 

01 

8.85456 

01 

•9,42396 

01 

9.9947F 

ni 

1.05496 0? 


1.10886 02 

1.1605E 02 

L.2Q90E 02 

L.2527E 

** n 

w it 

1.2904E 02 

1.3203E 

02 

1.3403E 

02 

1 .34596 

02 

1.345 IE 

02 

1.3452F 0? 

T 

3o0991E 02 

3.1050E 02 

3.1084E 02 

3.11106 

02 

3.1123E 02 

3.1133E 

02 

3.U43E 

02 

3.11446 

02 

3.11416 

02 

3.1137F C2 


3.1130E 02 

3.1123E 02 

3.11156 02 

3.11076 

02 

3.11006 02 

3-10946 

02 

3.10896 

02 

3.10876 

02 

3.10886 

02 

3.1092E 02 

HHQ 

B.9089E-01 

8.39196-91 

8.88236-01 

3 .874 96- 

-01 

8. 96486-01 

3.86&9t- 

-01 

a. 8655fc- 

■01 

a.. 365 i£ 

-01 

K.C659E- 

-01 

a.8673E-wl 


8. 86916-01 

8. 87126-01 

8.87356-01 

B.8757E- 

-01 

B. 97736-01 

8.8 7966- 

-01 

8.80 096- 

■01 

8.88146- 

-01 

8. 38136- 

-01 

8.3799E-01 

H 

3.L100E 05 

3.11596 05 

3.11936 05 

3.12196 

05 

3.12376 05 

3.124 7C 

05 

3. 17 52E 

C5 

3.1253F 

05 

3.12506 

05 

3.1246E .r 


3.12396 05 

3.12326 05 

3.12246 05 

3.12166 

05 

3.12086 05 

3.1202E 

OS 

3.11986 

05 

3.L196E 

05 

3. U96E 

05 

3,12016 05 

**aisr. 

ALONG X= 0.3700»***NEAP 

WALL RE=9.2340E 02 












z 

4.04966-04 

7.9573E-04 

1.25396-03 

L.9483E- 

-03 

2.91156-03 

4.1712E- 

■0 3 

5.75106-03 

7.67176-03 

9.951SE-03 

1.260BF-02 


1.56586-02 

1.9U8E-02 

2.2977E-02 

2,ri98E-02 

3.17706-02 

3.6715E- 

■02 

4.20666-02 

4.T872E- 

-0? 

5.42006- 

■02 

6.11036-02 

u 

4, 81166 01 

5.455CE 01 

5.9470E 01 

6.4B04E 

01 

7.02946 01 

7.58276 

01 

8.13556 

at 

8.68506 

01 

9.22926 

01 

9.7655F Cl 


1.02916 02 

1.08006 02 

1-12846 02 

1.17276 

02 

1.21206 02 

1.244 86 

02 

1.2692E 

02 

1.28246 

02 

1.7824E 

02 

1,28236 02 

T 

3-10396 02 

3.U05E 02 

3.11436 02 

3.11746 

02 

3. 11956 02 

3.1209E 

02 

3.12186 

02 

3.12216 

02 

3-1221E 

02 

3.1219E 02 


3.12146 02 

3-12086 02 

3.1200E 02 

3.1193E 

02 

3.1185E 02 

3.1179E 

02 

3. 1173E 

02 

3.1170E 

02 

3.1169F 

02 

3. 11756 02 

RHO 

a. 96336-01 

8.94436-01 

0,93346-01 

8,92476-01 

8.91856-01 

8.91456-01 

8.91216- 

01 

B, 91116-01 

8,91116- 

■01 

8.91186-01 


8.91326-01 

8.9L49E-0L 

8.91706-01 

8.91976-01 

0.92136-01 

8.92336- 

■01 

8,92486- 

01 

8.92596- 

-01 

8,92596-01 

8.9244F-01 

H 

3-11486 05 

3.12146 05 

3-125ZE 05 

3.12836 

05 

3-13046 05 

3.13196 

05 

3. U27E 

05 

3.1 33-OE 

05 

3.13316 

05 

3.132BF 05 


3.13236 05 

3.13176 05 

3.1310E 05 

3-13G2E 

05 

3.12956 05 

3-128 BE 

05 

3.1782E 

05 

3.1279E 

05 

3. 12796 

05 

3.12846 05 


*#DIST, ALONG X= 0.4200^'**»NEAR NALL Rt. ••.9y35E C2 


Z 


4.Z430E-0A 8c3274E-04 l,31L6E-03 2.0369E-03 3,0423F-03 4.3559E-03 6.0015C-Q3 7.9994S-i)3 U0367E-02 1.3122E-02 
1.62786-02 1.9851E-02 2.38296-02 2-81686-02 3.28566-02 3.7911E-02 4.3362E-02 4.925Or-02 5.56396-02 6.25976-02 


U 


4.446&E 01 5.04116 01 5.49896 01 5.9965E 01 6.5104E 01 7.0309E 01 7.55386 01 6.0765t 0* S.5970E 01 9.1132E 01 
9.6222E 01 1-0120E 02 1-0597E 02 1.1042E 02 1. 1444E 02 l,1792E 02 1.2071E 02 1.2256E 02 U2304F 02 1.2295E 02 


T 


3. 1074E 02 3.1146E 02 B.liSBE 02 3.1222E 02 3.12476 02 3.1264E 02 3.1275E 02 3,12816 02 3.12846 02 3.1283E 02 
3.1280E 02 3.127SE 02 3.1269E 02 3.12616 02 3.12546 02 3.L247E 02 3.1Z41E 02 3-12376 02 3.1235E 02 3.1241E 02 


54 


K.H0 

9, 0077 E-01 

8.9a69E-01 

8.9748E-01 

8.9649E-01 

8.9577E-01 

a,95Z8E-01 

B.9497E-01 

8.9480E 

-01 

a-9473E“01 

3.9476E-01 


B,9484E-0i 

8.9499E-01 

8-9517E-01 

8.9537E-QI 

8.95586-01 

8.9577E-01 

8. 9595E-01 

8.9608E 

-01 

8.9614C-01 

3.9597F-01 

H 

3.U83E 05 

3.L255E 05 

3.1297E 05 

3.1332E 05 

3. 1357E C5 

3.1374E 05 

3.13B5E 05 

3. 139 IF 

05 

3.1393E 05 

3.1392E 05 


3, 1389E 05 

3.1384E 05 

3.1378E 05 

3.1371E 05 

3.1364E 05 

3.1357E 05 

3.L351E 05 

3. 1346F 

05 

3.1344F 05 

3,1350F CF 

=('+DiST. 

ALONG X= 0.4700*+**NEAR 

WALL RE=8. 

7183E 02 








1 

4.4957E-C4 

6.8055E-04 

l.3846f-03 

2.1462E-03 

3.1993E-03 

4.5719E-03 

6.28 75E-03 

8,3659E' 

-03 

1.0324E-02 

1.3678F-02 


l,6940E-02 

2.0627E-02 

2.4721E-02 

2.9176E-02 

3. 3977E-02 

3,9139F-02 

4. 46 8 7 E- 02 

5 .06972-02 

5. 7104 P-02 

6,4l03F-r)2 

U 

4.0609E 01 

4.6162E 01 

5i0473E 01 

5.5188E 01 

6.C072E 01 

6.5032F 01 

7.0026E 01 

7.5035E 

01 

8.0041E 01 

8.5029F 01 


a,99?3E 01 

9.4841E 01 

9.9552E 01 

1.Q399E 02 

1.0B07E 02 

L.n69E 02 

1. 1473E 02 

U1696E 

02 

1.1810F 02 

1..1505E 0? 

T 

3.UQ8E 02 

3.1184E 02 

3.1229E 02 

3.1265E 02 

3. 1293E 02 

3.1313E 02 

3.1326F 02 

3. 1334E 

02 

3.1338F 02 

3.1339F 02 


3.1338E 02 

3.1334E 02 

3. 1329 E 02 

3.1322E 02 

3. 1316E 02 

3.13096 02 

3.1302E 02 

3.1297E 

02 

3. 129^^6 02 

3.1299E 02 

RHQ 

9,0460E-01 

9.0241E-01 

9.Q112E-01 

9.0006E-01 

8.9926E-01 

8.9370F^01 

8.9832 E-01 

8.9809E' 

-01 

8.9797F-01 

8.9794E-01 


8.9798E-01 

8.98Q9E-01 

8.9824E-01 

8.9842E-01 

8.9S61E-01 

8.988 lE-01 

8.9999E-01 

a .9915F' 

-01 

5.9925E-01 

8 .9909F-01 

H 

3.121TE 05 

3.1293E 05 

3,1338E 05 

3.1375E 05 

3.1403E 05 

3.1422E 05 

3. 14 36 E 05 

3 .1444E 

05 

3. 144SF 05 

3.1449E 05 


3.1447E 05 

3-1444E 05 

3.1430E 05 

3.1432E 05 

3.1425E 05 

3.141 QE 05 

3.1412E 05 

3.1407E 

05 

3.1403E 05 

3.1409E 05 


*=*DIST. ALONG X= 0. 5200****NEAR ViALL RE=8.49llE 02 . . 

Z 4.7024E-04 9.197AE-04 1.4452E-03 2.2387E-Q3 3.33S0E-03 4.7624E-03 6,54415-03 8.69945-03 1.1245E-02 1.4195E-02 

1.7562E-02 2.1359E-02 2.5B69E-02 3.0L41E-02 3.5056E-02 4.0323F-0? 4.S978E-02 '5.2033E-0? 5.3552E-32 6.55965-3: 

U 3.7700E Oi 4.2893E 01 4.6922E 01 S-1341E 01 5.5942E 01 6,D633E 01 6.5389E 01 7.0175E 01 7.4979E 01 7.9786F 01 

6.4577E 01 d.9321E 01 9.3947E 01 9.8348E 01 1.0244E 02 1.0616E 02 1-0937E 02 L.1192E 02 1.1357E 02 1.1390E 02 

T 3.U34E C2 3.1214E 02 3.1261E 02 3.1301E 02 3. 1331E 02 3.1352E 02 3.1367E 02 3.1377E 02 3. L382E C? 3.1385E 02 

3,l3a5E 02 3.1382E 02 3- 137SE 02 3.1373E 02 3.1367E 02 3.1360E 02 3.1354E 02 3.1348E 02 3..1344E 02 3.1348E 02 

RHQ 9,07a4E-01 9.0S52E-01 9.0415E-01 9.0301E-01 9.0Z16E-01 9.0154E-01 9.0111E-01 9,OQ83E-Ol 9.0066E-01 9.0059E-01 

9.0059E-01 g.Q066E-01 9.0078E-Q1 9.0094E-01 9.0112E-01 9.0130E-01 9.0149E-01 9.0165E-01 9.0178E-01 9.0167F-01 

H 3.1243E 05 3.1323E 05 3-1371E 05 3.1410E 05 3.1440E 05 3,1462E 05 3, 1477E 05 3.1487F 05 3. 1492E 05 3.1495E 05 

3.1495E 05 3.1492E 05 3. 1488E 05 3.1483E 05 3. 1476E 05 3.14T0E 05 3.1464E 05 3.1458E 05 3.1453E 05 3.1457E 05 




IS 


s 

& 


EX 

PRESS 


P-RECV 

WALL SHEAR 

HEAT FLX 

B.L. THK 

DISP THK 

MOM THK 

SHAPE 

FAC 

0.0 

7.100E 

04 

0,0 

3.967E 

01 

0.0 


1. 5536-02 

1,9136-03 

1,4716-03 

1.301E 

00 

3.906E-05 

7.100E 

04 

I.536E-04 

3.963E 

01 

7.372E 

00 

1.5546-02 

1. 9136-03 

1,4716-03 

1.301E 

00 

7.812E-05 

7.100E 

04 

2.583E-04 

3.96LE 

01 

1.37BE 

01 

1.554E-02 

1.9146-03 

1.471L— 03 

1.3016 

00 

l,562E-0<t 

7.101E 

04 

5.340E-04 

3.955E 

01 

2.7546 

01 

1.555E-02 

1.915E-03 

1.472E-:03 

1.301E 

00 

3. 125E“0<i- 

7.102E 

04 

1.0546-0 3. 

3.945E 

01 

5.365E 

01 

1.556E- 02 

1.917E-03 

1.473E-03 

1.30LE 

00 

6.250E-04 

7.104E 

04 

2.084E-03 

3.925E 

01 

1.0275 

02 

1.559E-02 

1.921E-03 

1.476E-03 

. 1.301E 

00 

1.250E-03 

7.L07E 

04 

4.X19E-03 

3,690E 

01 

1.8996 

02 

1.564E-02 

l,929E-03 

1,4816-03 

1.3026 

00 

2.500E-03 

7.X15E 

04 

a. X95E-03 

3,8296 

01 

3.364E 

02 

1.5746-02 

1.9456-03 

1.4926-03 

1.304E 

00 

5,UOOE-03 

7.129E 

04 

1.624E-02 

3.727E 

01 

5.637E 

02 

l,592E-a2 

X,9r8E-03 

1.513E-03 

. 1.308E 

00 

l.OOOE-02 

7.15fE 

04 

3. 189E-02 

3.565E 

01 

8.9066 

02 

1.625E-02 

2.0436-03 

L.555E-03 

1.3146 

00 

2.OO0E-02 

7.2UE 

04 

6.I76E-02 

3.306E 

01 

1.330E 

03 

L.690E-02 

2.1726-03 

1.6416-03 

1.3246 

00 

3.000E-02 

7.252E 

04 

a.508E-02 

3.159E 

01 

1.596E 

03 

1.8146-02 

2.2826-03 

1.7176-03 

1.3296 

00 

<t.Da0£-02 

7.294E 

04 

I.085E-01 

2.994E 

01 

1.3226 

03 

1.902E-02 

7.4016-03 

1.790E-03 

1.335E 

00 

5.000E-02 

7.333E 

04 

1.30XE-01 

2. 8626 

01 

1.997E 

03 

1,9706-02 

2.517E-03 

1,8786-03 

1.3406 

oo 

6.0006-02 

7. 37 IE 

04 . 

X.5X2E-01 

2.727E 

01 

2.1496 

03 

2.032E-02 

2.639E-03 

- 1.963E-03 

1-344E 

00 

7.0006-02 

7.405E 

04 

1.705E-01 

2.6166 

OL 

2.272E 

03 

2.092E-02 

2.759E-03 

2.047E-03 

1.3486 

00 

a, 0006-02 , 

7.439E 

04 

X.BS9E-0X 

2.5096 

01 

2.377E 

03 

2.1506-02 

2.8816-03 

Z.131E-03 

1.3526 

00 

9.000E-02 

7.470E 

04 

2.067E-01 

2.407E 

01 

2.46BE 

03, 

Z.265E-02 

3. 0076-03 

2.219E-03 

1*3556 

00 

i.OOOE-01 

7.50DE 

04 

2.234E-01 

2.313E 

01 

2-546E 

03 

2.3546-02 

3.1356-03 

2.3076-03 

1.3596 

00 

i.lOOE-OI 

7.523E -04 

2.391E-01 

2.229E 

01 

2.612E 

03 

2.427E-02 

3. 2626-03 

2.396E-03 

1.362E 

00 

1.2Q0E-01 

7.555E 

04 

2.542E.-01 

2. I486 

OL 

2.6&8E 

03 

2.493E-02 

3.3926-03 

2. 4866-03 

1.3656 

00 

1.300E-01 

7.581E 

04 

2.686E-0X 

2.Q71E 

01 

2.717E 

03 

2. 556 E- 02 

3.5246-03 

2.5776-03 

1.3676 

00 

1.400E-01 

7.6Q6E 

04 

2.823E-0X 

1.997E 

01 

2.7586 

03 

2.617E-02 

3.6606-03 

2.671E-03 

1.3706 

6 o 

1.5Q0E-O1 

7.629E 

04 

2.952E-0X 

1.931E 

01 

2.7926 

03 

2.678 6-02 

3,7946-03 

2.7646-03 

1.3736 

□ 0 

1.600E-01 

7.651E 

04 

3.076E-01 

1.867E 

01 

2.821E 

03 

2. 7766- 02 

3.9306-03 

2. 8586-03 

1.3756 

00 

i.700E-01 

7.672E 

04 

3.X96E-01 

1.805E 

01 

2. 8456 

03 

2,8786-02 

4.0716-03 

2.9546-03 

1.3786 

oo 

i.aooE-01 

7.693E 

04 

3.309E-01 

1.746L= 

01 

2.865E 

03 

2,9626-02 

4.213E-03 

3.0*526-03 

1.380E 

00 

L.yO06-01 

7.712E 

04 

3.4X6E-01 

1.693E 

01 

2.3816 

03 

3.034E-02 

4. 3546-03 

3. 1496-03 

1.3836 

00 

2.UO0E-OI 

7.73GE 

04 

3.520E-0X 

1.641E 

Q1 

2.393E 

03 

3.1016-02 

4.4976-03 

3.2476-03 

1.3856 

00 

2. lOOE-Ol 

7.748E 

04 

3.6X9E-01 

1.59 IE 

01 

2.902E 

03 

3.1666-02 

4,6446-03 

3,347E-'13 

l,’88F 

00 

2.200E-OX 

7.766E 

C4 

3.7X7E-0X 

1.5396 

01 

2.907E 

03 

3.23DE-02 

4. 796E-03 

3.449E-03 

1.3906 

00 

2.300E-01 

7.784E 

04 

3.322E-0X 

1.47ZE 

01 

2.904E 

03 

3.2956-02 

4.9696-03 

3.5626-03 

1,3956 

00 

2.400E-01 

7.801E 

•04 

3.-9X7E-0X 

1.421E 

01 

2.903E 

03 

3.3706-02 

5.1376-03 

3.6736-03 

1.3996 

00 

2.500E-01 

7.819E 

04 

4.014E-OX 

1.362E 

01 

2.8956 

03 

3*491£-02 

5.318E-Q3 
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APPENDIX A 


In the computation of confined flows the pressure gradient along 
the flow is usually the unknov/n parameter to be calculated. The finite 
difference equations given in section 2, along v/ith the supporting rela- 
tions given in Appendix A of Ref. 7, favor the case in which the pressure 
gradient is the specified varaible and the duct area along the flow is 
the unknown quantity to be calculated. When the pressure gradient is 
the unknown quantity these equations, as discussed earlier, are solved 
iteratively. Two questions arise in regard to the pressure iteration. 
One, how to adjust the previously used Ap (pressure drop across Ax) 
for the succeeding iteration. Second, how best to estimate starting Ap 
at the beginning of a pew integration step. In this appendix we have 
summarized our current know-how in regard to these tv;o quisstions. 

Applying conservation principles approximately across x and x + Ax 
we can derive the following equations: 


AiA2[i “ 


A‘ = -^P' 
m 


Equation (A-T) is derived in Appendix C of Ref. 7 and Eq. (A-2) follows 
easily from the equations given in that appendix. In these equations 


1)2 is average velocity at x + ax 
A-] duct area at x, 

A2 duct area at x+Ax 
iTi mass flow rate 
P2 pressure at x + Ax 
A’ small change in Ag 
p' small change in Ap 
U' omall change in Ug 
9 = infinity for incompressible flows 
= ratio of specific heats for compressible flows 

We also introduce the following nomenclature for later use: 

DPDTL specified tolerance within which the pressure gradient is 
to be calculated (put into the program as a dimensionless 
number, DPDTL = 0.01 means that the pressure graduent is to 
be calculated within +1%), 

ZTLRN tolerance for the calculation of the duct area, 

UTLRN tolerance in the calculation of velocity at x + Ax. 

At the beginning of the MAIN program DPDTL is specified through an 
assignment statement. Given DPDTL (p') Eq. (A-1) specifies the tolerance 
limit within which the duct area must be calculated to meet the pressure 
gradient tolerance. Equation (A-2) tells how closely Ug must be calculated 
during the velocity iteration loop. At the beginning of a new integration 
step Eqs. (A-1) and (A-2) are used to calculate ZTLRN (A') and UTLRN (U’) 
from the specified DPDTL (p'). Then with an initial estimate of Ap we 
start the velocity iteration loop until Ug at x+Ax converges within 


UTLRN. We next check calculated Ag against the specified Ag. If the 
difference, called A', is greater than ZTLRN we change Ap by p‘, 
calculated from Eq. (A-1) using A', for the succeeding pressure iteration. 
Note Eq. (A-1) does double duty— at the beginning of a new integration 
step it is used to calculate ZTLRN for that integration step and then 
during the pressure iteration loop it is used to correct Ap for the 
next iteration using A’ from the last Iteration. In the actual computa- 
tions Eqs. (A-1) and (A-2) are used with Ug replaced by the known center- 
line velocity at x, Pg by p^ , Ag by A^, and U' as the variation of Ug 
(U along streamline 2). To make up partly for these approximations and 
partly for those that went into the derivation of Eq. (A-1), Eq. (A-1) 
when employed within the pressure iteration loop, to calculate adjust- 
ment to Ap from A', is used along with a damping factor. Correction 
p‘ calculated from (A-1) is decreased by the factor 


DPDAHP = 0.5X 6XP (-0.115 ) • (A-3) 


We next turn our attention to the estimation of initial Ap, call 
it Apl, at the beginning of each new integration step. Denote the 
final calculated Ap at the given integration step by ApF. In our ear- 
lier computations Apl at x.^ was calculated from ApF at x^_g and x.-_i 
using straight line projection. In Fig. (A-la) is shown the distribution, 
using this approach, of Apl(x) and ApF(O) along part of a duct. The 
flow is through a rectangular diffuser. Note the flip-flop nature of 
the distribution of Apl with respect to mean Ap/Ax versus x curve. Note 


^ Initial estimates of Ap/AX 


©Final calculated values of Ap/AX 


X ^xa. 


Fjg. (A-1a). Initial pressure gradient estimated from 


previous two calculated values using a 
linear relation. 


Fig. (A“1b). Initial pressure gradient estimated from 


previous three calculated values using 
regression analysis. 
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also that to estimate Apl by fitting second order curve through ApF 
at ^‘i-l worsen the estimate of Apl. 

Since ApF is calculated not exactly but within a specified limit, 
the values of ApF will fall randomly about the mean Ap/Ax versus x curve. 
Thus Apl at x^. should be calculated from ApF at previous steps on stat- 
istical basis. In our later computations Apl at x-j was projected from 
ApF at x^_g, x^_ 2 » and x^‘_-j using regression analysis. Distribution of 
Apl and ApF using this approach is shown in Fig. (A-lb). Figures (A-la) 
and (A-lb) are for the same flow over the same section of the duct. Es- 
timation of Apl using regression analysis has improved the convergence 
rate in our computations. 

In the program listing given in section 3 Apl at x^. is calculated 
from previous three values of ApF using the following formulae: 

Ap = DPSLOP xi + DPZERO (A-4) 


DPSLOP = 


SSAp^X^ - (2Ap.)(EXi) 
3 - ( S 


(A-5) 


DPZERO = 


( £Ap^O(Sx.^) - ( Zxj)( SXjAp^-) 
3 2X;j^ - { Sx.)^ 


(A-6) 


In the program x at the first of the three points to be projected is 
taken equal to zero. Near the entrance before three ApF are available 
Apl is set simply equal to the last ApF. 
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ABSTRACT 

A computational method to calculate momentum and energy transport in two- 
dimensional viscous compressible duct flows is presented. The flow in the 
duct is partitioned into finite streams. The difference equations are then 
obtained by applying momentum and energy conservation principles directly to 
the individual streams. The method is applicable to laminar and turbulent 
flows. 


I. INTRODUCTION 


Two-dimensional effects in duct flows, such as the growth of momentum 
and thermal boundary layers, have been studied for a long time by separ- 
ating the flow into a core flow and a boundary layer [1]. Two-dimensional 
effects associated with sudden changes in duct shape in incompressible 
flows have been successfully computed using vorticity ti-ansport equation 
[2]. In high temperature compressible flows such as those encountered 
in magnetohydrodynamic (MHD) channels separation of flow into a uniform 
core and a boundary layer is difficult except at the very beginning of 
the channel. The difficulties arise from (1) large temperature differences 
between the wall and the core flow which lead to temperature gradients 
within the core and (2) shrinking size of the core as the flow proceeds 
down the channel. Also in calculating performance of diffusers with high 
inlet blockage factor, as in MHD diffusers, one cannot separate the flow 
into a core and a boundary layer. Computation of such flows have been 
recently carried out mostly using equations obtained by von Mises trans- 
formation of boundary layer equations [3,4,5]. Some MHO channel calcula- 
tions have also been carried out using Fluid-in-Cel 1 Method [6,7]. 

We have in this article presented a new approach to compute internal 
flows without separating the flow into a core and a boundary layer. A 
novel feature of this ’approach is the method by which the difference 
equations are obtained. In most of the approaches one starts with the 
differential equations, such as energy equation, Navier-Stokes equations, 
which are obtained form the conservation principles. These equations 
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either in their primitive forms or in transformed forms are then differ- 
enced to obtain the finite difference equations needed to compute the 
flow. In the present approach the flow in the duct is portioned into 
finite streams and the finite difference equations are then obtained 
by applying conservation principles directly to the individual streams. 
The resulting equations are similar to the ones obtained previously [3,4] 
by the procedure: derive boundary layer equations from the conserva- 

tion principles via energy and Navier-Stokes equations—carry out von 
Hises transformation of the boundary layer equations— finite difference 
the transformed equations using integral method finite differencing. In 
contrast, in the present approach there is direct relationship between 
the finite difference equations and their roots, the conservation 
principles. This direct relationship, we think, makes it easier to 
introduce approximations allovi?ed by the flow being simulated and also 
to incorporate any special features of the flow. 

In section II we have presented the finite difference equations 
based on the new approach and in section III discussed their solution 
techniques. The discussion is limited to two-dimensional flows. The 
application of the present approach to compute three-dimensional effects 
is straightforward and Is indicated at the end of section III. In the 
main body of the paper we have tried to present the new approach in 
broad and general terms to make its application easier to specific 
problems requiring individual handling. Some specific details are 
given in Appendices A and B. To provide an assessment of the present 
approach we have carried out sample computations of turbulent flow in 
a smooth pipe. The results are shown in Appendix C and are compared 
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with experitnental data and other works. We have also used this approach 
to compute two-dimensional viscous compressible turbulent flovis in 
hydrogen-oxygen combustion driven MHD ducts [8] and are currently 
using it to calculate diffuser performance with high inlet blockage 
factors . 
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II, FINITE DIFFERENCE EQUATIONS 

Consider two dimensional duct flow as shown in Fig. 1. The direction 
X is taken along the flow and z, measured from the duct wall, normal to it. 
For circular uucts z represents the radial distance from the wall. For 
rectangular ducts the variation of flow properties normal to z and x are 
neglected. Vie partition the total flow in the duct into a sequence of finite 
streams. To accomplish this partitioning we draw a series of streamlines 

1,2, as shown in Fig. 1. Streamline 1 is taken along the wall and 

streamline N represents the duct center or the symmetry line. The spacing 
between any two consecutive streamlines is determined by the desired reso- 
lution of the local gradients normal to the flow. We define streamline 
j + 1/2 as the streamline passing through the middle of j and j + 1. Stream- 
line j-1/2 is defined similarly. The stream j, in general, is now defined 
to be the flow bounded by streamlines S~ 1/2 and j + 1/2. The shape of the 
stream, of course, depends on the shape of the duct; in rectangular ducts 
the streams will be plane layers and in circular ducts the streams will be 
cylindrical shells. Stream 2 is defined to be the flow between the wall and 
streamline 2 + 1/2. Stream N is defined to be the flow between the stream- 
line N- 1/2 and the center or symmetry lifie N. 

We introduce the following nomenclature; 

ip. ; mass flow rate between the streamline j and j + 1/2. 

J 

; mass flow rate between the streamline j and j-1/2. 

J 

< 

Note + Tp. represents the mass flow rate through stream j. For calcula- 

J J 

tions of ip.'s see Appendix A. 

V 
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u,-(x), T^(x), p.(x), h.(x), C.(x); velocity, temperature, density, 

J J J J J 

enthalpy and constant pressure specific heat along the stream- 
line j. 


z.(x); distance of streamline j from the wall, 

J 

cross-sectional area of stream j (see Appendix A). 

J 




effective viscosity and thermal conductivity 


(including turbulent contribution when the flow is turbulent) 
along streamline j + 1/2. These are the coefficients that are 
later used to calculate momentum and energy transfer between the 
stream j and j + 1 across the interface at j + 1/2. 


Other nomenclature is introduced as we go along. 

The difference equations to calculate the flow properties are now ob- 
tained by applying conservation principles directly to the individual streams. 
We assume that the streamlines are chosen sufficiently close so that the varia- 
tion of the flow properties betv^een any two successive streamlines j and j + 1 
can be approximated as linear. Note the flow properties u., T., etc., though 
not their gradients, are continuous across these streamlines. However, across 
the interface of the streams, i.e. across streamlines j + 1/2, the gradients 
are also continuous thus determining, unambiguously, momentum and e^^ergy 
transfer between the adjoining streams. 

Applying momentum balance to stream j between x and x + AX, i.e. applying 
the momentum ba'^''nce to the control volume bounded by x, x + Ax, j-1/2, and 
j + 1/2, we obtain. 


Momentum Flux Out - Momentum Flux In 

at X + Ax 


at X 


= Pressure Forces + Viscous Forces + Body Forces 


0 ) 
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Estimation of the various quantities in Eq. (1) allows a great deal of flexi- 
bility to accommodate the desired accuracy and any special features of the 
given flow. We here present one such approximation for stream j, (streams 2 
and M require special treatment and arc considered later separately). 










^ v( 


3u • + u . 
J 


1 






On the left hand side of Eq. (2) the expressions involving ijj. represents 

V 

the momentum flux between j and j + 1/2 and those involving ti*." the momentum 

J 

flux between j - 1/2 and j. The superscript + indicates the unknown velocities 
at x+Ax, Ap the pressure difference between x + Ax and x, and F- the sum of 
body forces on the control volume. Also we have introduced 


''j+l/2 


‘ 


''j-1/2 


''.i-1/2 ,5 

Z3-Z..-/M/2 


(3) 


where ^/-l/2 upper and lower surface areas of the 

control volume (see Appendix A for calculation of A 's). Terms 
+ 

^j-1/2 defined similarly with y's and z's evaluated at x + Ax. Note we 
have estimated the viscous forces on a given surface by taking the average 
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of those at x and x + Ax. Rewritting i-q. (2) with all the unknown velocities 
collected on the left hand side we obtains 
+ 


+ r ’^i 1 H 


+ 

j+1/2 


[ 4'^j "'■ 4 ’^j 2'^j+l/2 2''j-1/2j 


+ r 1 + 

Vl 


fi 

4 ■ 2 'j-1/2- 


-apA. 


r 1 ]' 


^i 1 
4 ^2''j+l/2 


“j&f/ ■*■ fi'j' 


_ ly _ I 7 

2''j+l/2 2''j-l/2- 


“j-i ''j 


j = 3,4 N-1 


(4) 


Similarly, application of momentum balance to stream N, after some re- 
arrangement, yields 


1 , “ 4-lu + 

4% 2Vl/2 


- u 


N-1 


-I* 


fL+iw-^ 

■4 2^N-l/2 


= -ApA.^ + u, 


3 . - 1 




z'^n-y? 


+ u 


N-lL 


•I'l 


liL+lv 

4 2^i-l/2J 


+ F 


N (5) 


We note that Eq. (4) will yield Eq. (5) for j=N if we were to set 4'[y|^ , 

'^N+1/2’ ‘^N+l zero. Thus by initializing these quantities equal 

to zero at the beginning of a computer program stream N can be treated by the 
general formula for stream j. 

Next we will treat Stream 2. This stream requires considerably more com- 
putational finesse. The sources of the difficulty are (1) occurrence of steep 
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gradients nomal to the flow near the wall and (2) reckoning of wall surface 
roughness for turbulent flows. For laminar flows the wall surface roughness 
does not affect the flow and the gradients are not as steep and the flow 
can be simulated by taking streams sufficiently thin next to the wall. For 
turbulent flows we let go of the assumption that the velocity varies linearly 
between the wall and streamline 2. The momentum flux through the lower half 
of the stream 2 is estimated as some a times i{i 2 ~U 2 where a takes into account 
the nonlinearity of the velocity distribution next to the wall. Also the 
viscous force at the lower face of stream 2, that is at the wall, is calcu- 
lated from the Reynold's number along streamline 2 taking into account the 
wall surface roughness. Proceeding accordingly, we obtain by applying momen- 
tum balance to stream 2 



+ 1^2 otU2 


■] - 


SUj+Ug 


aU 


2J 


= -ApA2 -s* 2['^2-M/2^^3 ^2 ^ '^2+l/2^^3 ~ ^2^. 


2[Vi Ug + V^Ug] 


+ F, 


(6) 


Coefficients and Vg^-j^g are given by the general expression for 

presented earlier. For programming ease viscous forces along the wall have 

been expressed as V-jUg similar to viscous forces at the other interfaces. 
However is defined differently from other V's. Estimations of and a 
for turbulent flows are discussed in Appendix B. When the flow is laminar 

Eq. (6) can be used for the momentum balance for stream 2 with cx set equal 

to one half and V-j determined by 
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V 

1 n 


(Laminar Flow) (7) 


where the A-j^ represents the area of the lower surface of stream 2 and 
the viscosity, at the wall (V.^ is determined similarly). Rearranging 
Eq, (6) we obtain 


"*^3 [“ A ’*‘ 2 ^ 2 + 1 / 2- ^2 4’^2 “*^ 2 ’^ 2 '^ 2 + 1 / 2 "'^ 2 


*5" » 

= -ApAg^ + ^3 [“I "^^^2+1/2. “2-l’^2 ■*■ ‘^'*'2 


Equations (4), (5), and (8) provide N-1 equations to solve for the N-1 
unknown velocities Ug - 

X 

Equations for the unknown enthalpies h.‘ at x + Ax are obtained by apply- 

V 

Ing conservation of energy to the control volume under consideration; 


Energy Flux Out! - Energy Flux Inj 

•at X + Ax 'at X 


Heat Added. By Conduction + VJork Done by Viscous Forces + 

Q. (Internal Heat Generated) + W. (Work Done by the Body Forces) (9) 

J J 


The different terms in Eq. (9) can be estimated with varying degree of accur- 
acy (and the accompanying complexity) to reflect their relative importance 
for the flow being simulated. One estimate for stream j is as follows: 

Energy Flux In ^ -i- 3h^ + +^(3Uj^ + u ) 
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Expression for energy flux out is obtained from the expression for the 
energy flux in by adding superscript + to h's and u's. 


Heat added by conduction across the upper surface of the control volume 


where we have introduced 




*^ 3 + 1/2 




A? 


■j+1/2 


and 


M - J-^1/2 .s 

''j+1/2 " Cj{Zj+-,-Zj) ^j-M/2 


( 10 ) 


The H terms can be expressed in terms of V by using Pr, the Prandtl number. 
Thus we can v/rite 


H = ''j-H /2 where Pr - 


( 11 ) 


Work done by the viscous forces on the upper surface 


u ■ , T t u . 




) 




Expressions for the heat transferred and the work done across the lower 
surface can be written down from the corresponding one's for the upper 
surface by replacing j + 1 by j, j by j - 1, and j + 1/2 by j - 1/2. Substi- 
tuting the various estimates into Eq. (7) and rearranging to get all the 
unknown h^ 's on the left hand side of the equation we get. 
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r 1 

il — + + 

4 ^2^+l/2 
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for j = 3,4,. . . ,M-1 


( 12 ) 


The equation that results on applying the energy balance to stream N 
can be obtained from Eq. (12) by letting j equal to N and setting , 

*^n+l’ ^n+r '^N+r ^N+r ^N+1/2’ ’Wl/2’ ^N+1/2’ ''n+1/2 

Next, following the procedure outlined for the momentum balance, we 
apply energy balance to stream ? to obtain. 


1 1 




j ! 

i 


t 

j^CShg+hg) + |-(3 u 2 -i- U 3 )j 


h 
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^ ^2 ~ ^ ^ ^ 2 ^2 J ~ ^ 2^^2 ^ 

^ i [}*2+l/2^^3 " ^2 ^ ■*■ ^2+l/2^^3" ^2^] " ^w 

i [^2+l/2^^3^ " ^2^ ^ ^2+l/2^^3^“ “2 ^j '*' *^2 ^'^2 


(13) 


In Eq. (13) 3 and 6 are correction factors to take care of nonlinear varia- 
tion between wall and the streamline 2 when the flow is turbulent and are 
further discussed in Appendix B. Also in Eq. (13) h,, represents the enthalpy 
evaluated at the wall temperature and q the heat transfer to the wall. For 

W 

flows with wall temperature specified can be written as 

where is calculated from via the relation given in Eq. (11). Vihen q^^ 
is specified Eq. (14) is used to eliminate and appearing on the left 
hand side of Eq. (13). Equations (12), (13) and (14) provide H equations 
for the N unknowns hj^, hj^j_i , . . . jh^ and either one of or q ^ 

-J- ^ 

The proceeding equations for the calculations of u. 's and h. *s are to 

J J 

be supplemented with equations given ■'n Appendices A and B to calculate the 
duct width at x + Ax in terms of u. ‘s and p. 's. 

J w 
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In summary, Eqs. (4), (5). (8). (12), (13), (14), (B-2), (A-13) 

[or (B-3), (A-14), and (A-15)], along with the thermodynamic relations 

i 

interrelating hj p, T, and p are our basic set of equations to calculate 
the two dimensional duct flow. 


Ill, SOLUTION TECHNIQUE 


Overall methodology adopted to solve the duct flow equations will 

depend upon the given design specifications, e.g. pressure gradient is 

specified or variation of duct width is specified, etc. In any case one has 

to solve the momentum and energy difference equations given in section II. 

'*{* 

The set of simultaneous equations to solve for new u. , as also the set 

J 

+ 

to solve for new hj , are of the general form 

-a.X.,1 + b.Xj - = dj. (15) 

where X. stands for the unknowns u. or h. . Coefficients a., b. and c. 

J. 0 • J 3 J J 

+ + 

due to their dependence on V. or H. are indirectly functions of X.. 

J J J 

Thus Eq. (15) is nonlinear. If we had approximated the viscous forces 
and the heat transfer across stream interfaces by their values at x rather 
.than by their averages between x and x + Ax coefficients a-, b. and c* would 

cl J l1 

have been independent of X.. In that case Eq. (15) will constitute a set 
of linear tri-diagonal equations and can be solved directly by the standard 
algorithm to solve such equations. In the present case these equations are 
solved by iteration. 

At first we calculate a., b., and c. for the momentum difference equa- 

Jo J 

tions setting equal to V and solve for new using the algorithm to 

J 

+ + 

solve directly tri-diagonal equations. With these new u. and H set equal 

J 

to H we calculated a., b., c- and d. for the energy difference equations 

J J J J 

and solve for h. . Next we calculate z- using newly calculated flow prop- 

J J 

erties at x + ax. The procedure is then repeated this time around calculating 

"f" ’J" 4* 4* 

V and H using the newly calculated values of u- , z. , etc. The iteration 

J J 
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1s continued until the desired convergence is achieved. The number of required 
iterations increases^ in regions of rapidly changing temperature and velocity 
profiles. These are also the regions where had we calculated a-, b. and c. 
from V and H alone (instead V, V and H, H ) the error introduced by such 
approximation would have been the greatest. 

When the pressure gradient is specified the z^^'s calculated at the end 
of the last iteration give the channel width at x + ax. When the channel 
vddth is specified iteration on pressure gradient. is required. The itera- 
tion in the pressure gradient can be performed either over and above the 
Iteration on V and H or combined into the V -H iteration loop. In the 
latter case at the end of each iteration v;hen new z. and thus the channel 

J 

width at x+Ax is calculated the assumed pressure gradient is readjusted. 

We conclude this section by indicating how the present approach can be 
used to calculate three-dimensional flows. In three-dimensional case one 
win partition the flow in the duct into three-dimensional streams. In rec- 
tangular channels the streams now instead of being plane layers as in tv/o- 
dimensional case will be rectangular bars. The difference equations to 
calculate the flow can now be written down as in the two dimensional ca 
except now when calculating viscous forces acting upon the control volume 
and the heat conducted into it one vnll have to take into account four neigh- 
boring streams. The resulting equations now will be pentadiagonal and will 
require more computational time to solve. 
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APPENDIX A 


In this appendix we present formulae for calculating lij-”, A.^, A-^, 

J J J J 

and z. for the rectangular and cylindrical ducts. These are the formulae 
we have used in our computational works and can be easily modified to fit 
the requirements of other problems. For cylindrical ducts we Introduce r- , 

J 

the distance of streamline j from the centerline. For duct of radius Rg, 
r. and z. are related by 

J J 


'■j ‘ "o - 


(A-1) 


For rectangular ducts the formulae are for a duct of unit v/idth. (Formulae 
are presented in forms to reveal their origin. They can be further contracted 
or simplified for computational efficiency.) 

MASS FLOW RATES 
(a) Rectangular Ducts 




J = 2,3,. .. ,W-1 


(A-2) 


- . / SPi-UjH-Pj-iUj-i 

I 4 2 




j = 3»4,... ,N 


(A-3) 
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(b) Cylindrical Ducts 




j - 2,3,. . . ,N-1 

(A-4) 



j = 3,4,.. .,N-1 

(A-5) 

73pf^u^+P,^„lUj^_-|\ 2 / 

% 4 

(A-6) 

Calculations of are discussed in Appendix B. 


CROSS-SECTIONAL AREAS 


(a) Rectangular Ducts 


A2 2 

(A-7) 


(A-8) 


J 2 

j = 3,4,...,N“1 


X _ " ^N-1 

‘N ' 2 


(A-9) 
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(b) Cylindrical Ducts 



(A-lOl 


(A-11) 

3 ~ 3,4, .. . IN’-! 


'■n-1 /‘^ 

(A-12) 


z. AND r- TERMS 
3 J 

When the p's and u's have been calculated at a new position along the 
duct the z's and r's at the new position are then calculated by the follov/- 
ing formulae. 


(a) Rectangular Ducts 

After Z 2 has been calculated as discussed In Appendix B we calculate 

Z 3 , 




= "j-l + 


(A-13) 


(b) Circular Ducts 


= b 




N-1 L ir(3pfjUf^+ ^PN-l“N 


:?r. 


1/2 


. (A-14) 


18 




+ 


4 ^,- 


|l/2 


Foa: ^ + uT ■*' "Ti:3u, +V.,)' '■j 


J - 


Radius is calculated as discuss-d in Appendix B. 


SURFACE AREAS 

(a) Rectangular Ducts 


for all interfaces 



(b) Circular Duct 


for all interfaces ^^+ 1/2 “ ^ ■** '^^>1 ) 


(A-15) 


(A-16) 


(A-17) 


APPENDIX B 


f 

In this appendix we will discuss how to calculate the various correction 
factors that take into'account the non-linear nature of velocity distribu- 
tion between streamline 2 and the wall. There are many approaches possible 
to calculate these correction factors. The discussion that follows presents 
a method rather than the only one. For flow through circular pipes thick- 
ness of stream 2, i.e. 22 -z-], is assumed so small compared with the duct 
radius that we can neglect the curvature of stream 2. In evaluating the 
various integrals we have held p constant equal to some p”- Fof incompres- 
sible flows F is, of course, just equal to p. For compressible flov/s p 
is taken- to be the average of p between streamline 2 and the wall. 

We introduce mass flux correction factor, y, by 

r-^J^ «dz (B-1) 

In terms of y we can write 

ipf =yFugZ2 (B-2) 

Equation (B-2) is used at the beginning to calculate Tp 2 ~‘ and later on 
to calculate Z 2 for rectangular ducts. For circular ducts r-j i> calculated 
by 

r-] = rg + V/'*'^2^2 


From their deifinitions we have for the correction factors a and 6, 


(B-4) 


?2 


a T^2 ^2 


- I 


p dz , 


that is 
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a = — I — / dz , 

■ 0 

and similariy, 

6 = — f u^dz (B-5) 

YU2 22 -^0 


If one assumes as 1s often done, that the variation of h in the near wall 
region is similar to u then the correction factor 3 will be obviously equal 
to oc. If some other variation of h is assumed then the integral 


3 = 





evaluated with the assumed variation of h will yield g. 

The velocity distribution in the near v/all region, needed to evaluate the 
proceeding integrals, ■’’a calculated from the relation 


du „ 

'^tdz ■ ■'w 


(B-6) 


where is the shear stress at the wall and is the viscosity including 
the turbulent contributions. To proceed further one must select some ex- 
pression for To illustrate the procedure we use von Driest 's relation [8] 
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for vt (this relation has been successfully used in our own computations [7] 
and by several othen authors [3,4]). For smooth walls we have from [8] 


11 ^ = y + pk^z^ ^ - exp(-z/f^/yA) ^ 


(B-7) 


where y is the laminar viscosity and k and A are constants. Introducing 
non-dimensional variables 


= z/fiijjjp/y , and u^ ~ pu/ 


(B-8) 


we obtain on integrating (B-6) using (B-7) 


r 

{ 1 + f 1 + 4k2z* 


- exp(-z*/A))2j’ 


(B-9) 


In terms of the non-dimensional variables integrals (B-1), (B-4), and 
(B-5) become, omitting the subscript 2, 


-^1 


— ^ f u/dz* 

TU* z* Q 


— I — / 

y^-k 2 * 0 
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Introducting friction factor, f, and Reynolds number. Re, we have in 
terms of non-dimensional variables 

I 

f = \ , and (B-13) 

Re = = u*z* . (B-14) 

After having calculated the velocity distribution from (B-9) y, a, 6, 
and f can be next evaluated as functions of Re once and for all and the 
results then later used during the flow field computations. We have numer- 
ically integrated these integrals us-ing k=0,4 and A = 26, the values re- 
commended by Van Driest [8] for smooth walls. The results of the cal- 
culations are shown in Figs- (B-1) and {B-2). The numerically obtained 
relation between any of the factors and Re can now be stored for later use 
'either by fitting a high order polynomial over the whole range of Re of 
interest or by breaking the range of Re into small segments and then ap- 
proximating the relation over each segment by a straight line. The quant- 
ity V.J introduced in section II in Eq. (6) is in terms of f given by 

Pi Po 

= f g— U2 (B-15) 

and V-j is given similarly with the quantities in the right hand side of 
(B-15) evaluated at x + AX- 
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APPENDIX C 


In this appendix we present two sets of computations of turbulent 
flow in smooth pipes using the present approach. In these computations 
shear stress at the v;all was calculated from the local Reynolds' number 
as discussed in Appendix B. Turbulent shear stress, pts in the main 
body of the flow was calculated using Prandtl's mixing length hypothesis. 

Variation of the mean free path, along' 2 , the radial distance from 
the v/all , was calculated from the equation 

I = 0.14- 0.08^1 0.06^1 (C-2) 

where R is the radius of the pipe. Relation given in (C-2) is obtained 
from the experimental data for pipes as discussed in Ref. [1], 

DEVELOPING FLOW 

In the first set of the computations we calculated the developing 
velocity profiles in the inlet region of a, smooth pipe starting with 
slug flow (radially uniform velocity) condition at the inlet. Compu- 
tations were carrird out with twenty radial grid points, with the grid 
spacing increasing as we moved away from the wall. Axial integration 
step, Ax, was taken at start equal to R/256. This Ax was doubled after 
each integration step until it reached the value R/4; from there on 
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Ax was held constant. The iteration in the pressure gradient was carried 

out over and above the iteration on V {refer to section III). The for- 

mula used to iterate on pressure gradient is discussed at the end of 

+ 

this appendix. The V iteration was carried until the velocity next 
to the vmll varied less than 0:02 percent, and the pressure iteration 
was carried out until the pipe diameter at the new x was within 0.02 

precent of its specified constant value. For each integration step 

- + 
it took a combined total of 5 to 6 Iterations to meet both V and pipe 

diameter convergence criteria. The initial guess for pressure gradient 

at each new x was obtained from the previous two values using a linear 

relation. 

The results of the computations are shown in Figs. (C-l) and (C-2). 
In these figures ut, represents the bulk velocity, p the density, r the 
distance from the center line, x the distance from the entrance, and 
D the pipe diameter. In Fig. '(C-1 ), we have plotted the computed radial 
vellcity profiles at three values of x/D. Also shown in this figure 
are experimental values obtained by Barbin and Jones [10]. Small differ- 
ences betv/een the measured and the computed values are, we think, due to 
the following reason: in the computations we used for the radial varia- 
tion of s, Eq. (C-2), which is for fully developed turbulent flow; in 
the experiment, as noted by the authors from their turbulent intensity 
measurements, in the inlet region the turbulent shear stresses near the 
center of the pipe are below their fully developed values whereas near 
the v/all they exceed their fully developed values. In Fig. (C-2) 
we have shown the computed pressure drop in the entry region. Also 
shown In this figure are the experimental values obtained in Ref. [10]. 


25 


The computed pressure gradient after about fifteen pipe diameters becomes 
essentially constant and is equal to the pressure gradient calculated 
using Moody chart [1]. 

FULLY DEVELOPED VELOCITY PROFILES 

In the second set of computations we calculated velocity profiles 
for fully developed flow in a smooth pipe. These profiles ware obtained 
by setting momentum flux out of the control volume equal to the momentum 
flux into the control volume. This can be accomplished easily, although 
not most economically with respect to compute ting time, by setting all 
the iji.'s equal to zero in the program designed to compute developing 

J 

flows. The results of the computations are shown in Fig. (C-3). The 
computed values are shown by discrete points to indicate the logarithmic 
nature of the radial grid spacing using in these computations. The non- 
dimensional parameters used in the figure are defined by 


z . 

V 


and 



pu 


where is the shear stress at the wall. Also shown in, this figure 
are van Driest's calculations and the relation u* = 5.5 + 2.5 s,u z*, both 
of which compare well with the experimental data. 

PRESSURE GRADIENT CORRECTION FORMULA 


As discussed in section III in the present approach when the pressure 
gradient along the duct Is not specified it is obtained by iteration. 

The area of the duct at x + Ax is calculated during the computations 
using a guessed Ap, the pressure drop across ax. The guessed Ap is then 
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iterated upon until the calculated area at x+Ax meets the requirements. 

To better guess Ap at the beginning of each new Iteration a formula, 
based on approximate application of conservation principles, has been 
derived. We denote area, average density, and average velocity by 
A-j , Pp at X, and by Ag, Pgj Ug at x+Ax. Applying mass and momentum 
balance across ax ’ we obtain 

» 

A^PlU^ = A^PgU' (C-3a) 

= -Ap-Ai-F^ (C-3b) 

v;here m represents mass flow rate and the viscous forces. Now we ask 
the question that to produce a small change in Ag how much change In Ap'* 
is required? Let Ag, p 2 > Ug, be another set of values at x+Ax that 
also satifies Eqs. (C-3) and differ from Ag, Ug, pg, Pg by only small 
amounts. We set 

Ag^Ag + A', = Ug + U% Pg = pg + p", p^ = pg + p" (C-4) 

Substituting relations from (C-4) into (C-3) we obtain after neglecting 
second and higher order terms in small quantities A", p"*, etc.. 



In obtaining Eq, (C-5) we have assumed (a) the small changes in Ag, pg, etc., 
induce negligible change in Fy, and (b) changes p’' and p" are related 
by 
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£1 = e (C^6) 

P2 P2 

where e is some constant. Equation (C-5) tells, approximately, that Ap 
when incremented by p-* will produce a change equal to A-' in Ag. During 
computations at the end of each pressure iteration the difference be- 
tween the computed and the specified area is calculated, calling it A" 

Eq. (C-5) is solved for p’'. The p"*, thus obtained, is in turn used to 
adjust Ap for the next iteration. In our experience so far, the most 
satisfactory value of e seemed to depend on the flow being computed. 

In the sample calculations given in this appendix, of incompressible 
flow, 6 was set equal to infinity. 

The sample calculations were carried out using FORTRAN IV programming 
language on IBM 360/44 with DOS operating system. For the case of devel- 
oping flow, discussed in this appendix, the program required 34 K of 
core memory. The computational CPU time was 2 minutes and 22 seconds to 
carry out the calculations to 60 pipe radii (approximately 240 integration 
steps). 
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